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SUMMARY 


Criteria  are  stated  for  the  existence  of  solitons  in  one-dimensional 
organic  conductors,  particularly  trans  polyacetylene.  All  physical  situations 
considered  are  described  in  terms  of  a  classical  anharmonic  oscillator  moving 
in  an  anharmonic  double  potential  well.  As  the  soli ton  is  known  to  exist  in 
the  line  of  pendula  joined  with  elastic  coupling  and  also  in  the  long  Josephson 
junction  in  the  form  of  the  fluxon,  the  theoretical  position  is  taken  that  if, 
indeed,  the  soliton  exists  in  polymers,  whether  of  the  vibrational  or  the 
charge-domain-wall- type,  then  it  should  conform  to  certain  established  criteria 
of  behavior  already  rigorously  studied  in  the  aforementioned  situations. 


In  order  to  understand  the  dynamics  of  a  conducting  one- dimensional  sys¬ 
tem,  whether  or  charge  or  mechanical  energy,  a  microscopic  parametric  excitation 
model  formalism  must  be  adopted.  On  the  other  hand,  in  order  for  this  model  to 
be  related  to  physical  measurements  such  as  that  of  current  and  voltage,  a 
microscopic  formalism  (the  Hubbard  model)  must  also  be  achieved.  Both  pictures 
have  advantages  and  disadvantages.  However,  both  pictures  are  necessary  for 
complete  understanding. 
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SOLITON  PROPAGATION  IN  ONE-DIMENSIONAL  ORGANIC  CONDUCTORS 

T.W.  Barrett 

Naval  Air  Systems  Command 
Code  AIR-031 
Washington,  D.C.  20361 

I.  PBEMHTNABTFR 

The  notion  of  a  propagating  soliton  is  not  a  simple  notion  and  the 
requirements  for  its  existence  are  rigorous.  The  aim  of  this  paper  is  to 
indicate  the  complexity  involved  in  its  existence  and  to  suggest  criteria 
for  such  existence.  Another  aim  is  to  indicate  the  fundamental  importance 
of  the  concept  of  soliton,  first  called  such  in  1965  (Zabushy  and  Kruskal, 
1965;  Zabusky,  1981),  but  whose  origins  stem  from  observations  by  John  Scott 
Russell  (1834).  The  timeliness  of  this  endeavor  is  due  to  the  consideration 
being  given,  on  the  one  hand;  to  proteins  (Davydov,  1977-1979;  Scott,  1981; 
Hyman  et  al.,  1981)  and  deoxyribonucleic  acid,  DNA,  (Balanovski  and 
Beaeonsfield,  1982)  as  mediators  of  vibrational- type  solitons,  and,  on  the 
other  hand,  to  trans  polyacetylene,  (CH)z,  as  a  mediator  of  soliton 
dependent  electrical  conductivity  or  charge-domain-wall  movement  (Su, 
Schrieffer  and  Heeger,  1979,  1980). 

The  charge-domais-wall-type  soliton  of  (CH)^  has  received  much 
attention.  The  case  for  the  existence  of  such  solitons  in  (CH)z  has  been 
made  mainly  on  the  basis  of  electron  spin  resonance  (ESR)  data  (Heeger  and 
MacDiarmid,  1981).  That  solitons  do,  in  fact,  exist  in  trans  polyacetylene 
could  be  treated  as  a  hypothesis.  On  the  other  hand,  the  mathematical 
entity  of  the  soliton  has  undergone  considerable  evolution  and  its  history 
has  been  described  by  Scott  (1981).  The  intention,  here,  is  not  to  question 
whether  the  mathematical  entity  of  soliton  propagation  can  exist  in  polymers 
(undoubtedly  it  does)  but  to  refine  judgment  concerning  when  it  can  and 
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when  it  cannot.  This  is  s  matter  of  discrimination  based  on  both  physical 
and  mathematical  reasoning. 

All  physical  situations  considered  here  are  described  in  terms  of  a 
classical  anharmonic  oscillator  moving  in  an  anharmonic  double  potential 
well  (DPW) .  As  the  soliton  is  known  to  exist  in  the  line  of  pendula  joined 
with  elastic  coupling  and  also  in  the  long  Josephson  junction  as  the 
fluxon,  the  attitude  has  been  taken  here  that  if,  indeed,  the  soliton  exists 
in  polymers,  whether  of  the  vibrational  or  charge-domain-wall-type,  then  it 
should  conform  to  certain  established  criteria  of  behavior  already 
rigorously  studied  in  the  aforementioned  situations. 

The  main  theme,  here,  is  that  in  order  to  understand  the  dynamics  of  a 
conducting  one-dimensional  system,  whether  of  charge  or  mechanical  energy,  a 
microscopic  parametric  excitation  model  formalism  must  be  first  achieved 
(seotions  I-IX) .  On  the  other  hand,  in  order  for  this  model  to  be  related 
to  physical  measurements  such  as  that  of  current  and  voltage,  a  macroscopic 
formalism  (using  the  Hubbard  model)  must  also  be  achieved  (sections  X-XI). 
Both  pictures  have  advantages  and  disadvantages.  However,  both  pictures  are 
necessary  for  a  complete  understanding. 

In  section  III,  the  possibility  of  sending  solitons  through  joined  one¬ 
dimensional  chains  into  one  arm  or  both  arms  of  a  T-junction  is  examined. 

The  conclusion  of  this  section  is:  if  it  were  possible  to  synthesize  a 
polymer  made  up  of  joined  segments  (of  CH)^  for  example)  the  synthesized 
polymer  would  demonstrate  different  phonon  modes  in  the  vicinity  of  the 
junction,  due  to  different  phonon-dependent  propagation  energies.  The 
likelihood  is  that  the  soliton  would  not  progress  past  the  junction  into  the 
arms  and  probably  would  not  approach  within  the  vicinity  of  the  junction. 

In  section  IV  the  concept  of  tunneling  is  related  to  a  parametric 
excitation  dynamic.  The  advantage  in  this  is  the  subsumption  of  tunneling 
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behavior  under  energy  conservation  rules,  rather  than  probabilistic  lavs. 
This  aim  is  of  importance,  because  the  mathematical  entity  of  the  soliton  is 
itself  rigorously  restrained  by  conservation  of  energy  lavs  (Bullough  and 
Caudrey,  1980). 

In  section  V  the  effect  of  pinning  solitons  is  addressed.  This  amounts 
to  the  trapping  of  a  discommensuration  so  that  it  is  no  longer  moving. 

Thus  by  pinning,  a  soliton  ceases  to  exist  (because  it  is  a  dynamic  and 
moving  vave  by  definition)  and  the  resulting  static  discommensuration,  or 
kink,  takes  its  place.  The  research  literature  has  sometimes  not 
distinguished  betveen  a  soliton,  vhich  is  a  moving  vave,  and  a  trapped  or 
pinned  discommensuration,  vhich  is  a  kink. 

In  section  VI  the  energy  conservation  lavs  for  the  order  parameter, 
energy  and  damping  constant  derived  for  nonlinear  systems,  are  used  to 
characterize  the  four  forms  Ohm's  and  Kirchoff's  lavs  vould  take  in  such 
nonlinear  or  soliton  exhibiting  systems. 

In  section  VII  the  dynamic  form  factor  or  spectral  density  is 
introduced  vhich,  for  a  system  in  a  displacive  transition  mode  such  as  a 
soliton  exhibiting  system,  contains  a  soft  phonon  vhich  pumps  the 
displacement.  It  may  be  emphasized  that  such  a  concept  is  out  of 
consideration  in  a  static  or  even  a  many  body  approximation  approach  to  one¬ 
dimensional  conduction.  The  dynamic  form  factor  does  not  even  appear  in  the 
macroscopic  landscape.  In  treating  the  dynamic  form  factor  the  advantage  of 
a  monomer  or  microscopic  treatment  of  solitons  is  thus  demonstrated. 

Hovever,  in  section  XI  and  XII  this  microscopic  treatment  is  relinquished 
for  a  macroscopic  or  Hubbard  many  body  approximation  so  that  long  chain 
systems  may  be  treated.  The  dynamic  form  factor's  soft  phonon  mode  remains 
a  facet  of  soliton  existence  criteria  and  of  the  microscopic  approach. 
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In  section  VIII  the  necessary  and  sufficient  conditions  for  soliton 
propagation  are  addressed.  Section  IV  contains  a  description  of  soliton 
propagation  in  a  system  of  changing  and  incomplete  degeneracy.  A  system 
capable  of  demonstrating  soliton  propagation  is  shown  in  section  X  to  be  one 
in  which  the  equation  of  motion  for  the  system  exhibits  broken  symmetry. 

This  relation  is  described  in  terms  of  renormalization  group  theory. 

Finally,  in  section  XI  a  macroscopic  picture  (the  Hubbard  model)  is 
introduced  and  related  to  the  double  potential  well  (DPW)  enharmonic 
oscillator  model. 

II.  INTRODUCTION 

It  is  well  known  that  the  interplay  of  currents  and  voltages  in 
Josephson  circuitry  is  similar  in  form  to  that  encountered  in  certain 
mechanical  systems,  and,  in  particular,  to  the  simple  pendulum.  The  ability 
to  understand  and  visualize  the  behavior  of  Josephson  junctions  in  terms  of 
mechanical  devices  has  greatly  facilitated  understanding  (Fulton  1975, 

1977).  As  (a)  the  fluxon  of  the  long  Josephson  junction  is  a  soliton 
propagation,  and  (b)  the  soli ton- with- spin  and  the  soliton-with-charge  have 
been  implicated  in  conduction  experiments  on  undoped  and  doped  trans¬ 
polyacetylene,  the  methods  used  to  study  the  Josephson  junction  fluxon  are 
readily  applicable  to  the  study  of  the  polyacetylene  soliton,  if  it  exists. 
The  suggestion  that  vibrational  solitons  exist  in  proteins  and  DNA  may  also 
be  studied  from  this  theoretical  platform. 

The  models  constructed  for  the  study  of  the  Josephson  effect  include 
those  of  Scott  (1969),  Valdram  et  al.  (1970),  Sullivan  and  Zimmerman  (1971), 
Yamashita  et  al.  (1974)  and  Fulton  (1977).  These  models  have  been  of  great 
use  in  developing  physical  intuition  concerning  the  effect  and  agree  in 
their  choice  of  mathematical  and  physical  description  of  that  effect.  This 
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agreement  means  that  actual  physical  model  bnilding  is  no  longer  necessary 
for  predicting  new  behavior  and  analysis  may  proceed  by  numerical  methods 
based  on  those  models. 

The  advantages  of  numerical  analysis  based  on  mechanical  models  are 
many.  Firstly,  the  performance  of  the  model  can  be  readily  altered  by  a 
change  of  variables.-  By  dimensional  analysis  polymer  compounds  for  possibl 
synthesis  may  be  suggested.  This  approach  is  more  systematic  than  the 
trial-and-error  method  of  obtaining  compounds  by  difficult  synthesis 
procedures,  testing  them  and  then  attempting  to  explain  their  function. 

Secondly,  recent  theoretical  approaches  to  the  polyacetylene  soliton 
emphasize  the  static  aspects  of  the  problem,  obscuring  the  dynamic  and 
diffusive  aspects.  An  analysis  based  on  a  physical  explanation  refocuses 
attention  on  these  neglected  aspects  of  the  problem. 

Thirdly,  there  are  a  number  of  variables  of  soliton  propagation 
indicated  by  mechanical  analysis  which  are,  as  yet,  unstudied. 

Hie  known  and  agreed  upon  conditions  for  fluxon  (soliton)  propagation 
are  now  examined  in  the  long  Josephson  junction,  and  these  are  used  to 
stipulate  the  proper  conditions  for  the  appearance  of  solitons  in  polymers. 
The  rationale  for  this  theoretical  posture  is  that  the  fluxon  (soliton)  in 
the  long  Josephson  junction  is  firmly  established,  but  the  polymer  soliton 
is  presently  a  theoretical  inference.  The  soliton  being  a  mathematical 
entity,  at  a  level  of  abstraction,  the  conditions  pertaining  to  the  former 
soliton  (the  fluxon)  must  also  apply  to  the  latter  (the  soliton  in  polymers 
whether  vibrational  or  moving  charge  wall).  Examination  of  the  soliton  of 
the  simple  pendulum  line  is  also  instructive  in  that  its  behavior  is 
described  in  familiar  physical  variables.  Thus,  the  intention  is  to  prop 
the  theoretical  conditions  for  polymer  soliton  propagation  upon  the 
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established  conditions  for  soliton  propagation  in  the  long  Josephson 
junction  and  the  pendnlnm  line. 

The  phenomenological  equations  describing  the  macroscopic  behavior  of 
superconducting  barriers  are  (Josephson,  1965): 

3<J>/3x  -  (2ed/ic)H  ;  3<|>/3t  *  (2 e/£)Vi  J  »  j,sin4>  +  cV,  (1) 

where  <J>  is  the  phase  difference  between  two  superconductors,  H  is  the 
magnetic  field  in  the  barrier,  V  is  the  voltage  across  the  barrier,  d  =  2\  + 
1 ,  A  is  the  penetration  depth  and  l  is  the  barrier  Sickness.  If  the 
barrier  is  regarded  as  having  a  capacitance  C  per  lit  area  and  using 
Maxwell's  equations  and  Equ.s  (1),  above,  we  obtf 


c*  at*  e* 


where  ca  *  ca/4ndC  is  the  phase  velocity  in  the  barrier.  A*  =  hc*/8jredj^ 
is  the  penetration  depth,  and  0  «  4ndcaa/ca  *  a/C  is  a  damping  constant. 


If  the  barrier  is  assumed  to  be  of  unit  length  in  the  y  direction,  then 


the  total  free  energy  of  the  barrier  is  (Lebwhol  and  Stephen,  1967): 

F  f dx [ ( 1— cost)) )  1  +  1/2  Aa  <|V  +  1/2  <^)*  (|$)*1,  (3) 

it  O  dx  C  dt 


where  the  first  term  is  the  coupling  energy  for  two  superconductors  and  the 
others  the  electromagnetic  energy  in  the  barrier.  The  rate  of  dissipation 
of  energy  is: 

dF/dt  *  -o/Vadx  (4) 
We  may  note  here,  in  anticipation  of  conclusions  to  follow,  that  removal  of 
any  terms  in  the  left  hand  side  of  Equ.  (3)  will  reduce  energy  dissipation. 

Recognition  of  an  analogy  came  about  when  Young  (1964)  pointed  out  that 
solutions  of  Equ.  (2)  representing  vortex  lines  in  the  barrier  are  obtained 
as  solutions  of: 

aa$/9xa  -  (1/^*)  »ia$.  (5) 
which,  except  for  sign,  is  the  equation  of  a  pendulum. 
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Next,  the  second  model,  the  pendulum  line,  is  considered  as 
mechanical  analog  of  the  long  Josephson  junction.  Figure  1  sketches  a 


simple  model  with  moment  arm,  1^,  moment  of  inertia,  M,  provided  by  a  fly 
wheel,  damping  coefficient,  D, 


Ic  sin* 


Fig.  1  After  Fulton  (1977). 

provided  by  a  paddle  wheel  arrangement,  and  phase  angle  <!>  measured  from 


vertical  downward.  The  rod  is  subject  to  an  external  torque  1. 

In  the  terms  of  such  a  model,  one  may  list  major  corresponding 
quantities  for  Josephson  circuits,  the  pendulum  analog  (Fulton,  1977)  and 
trans  polyacetylene : 

TABLE  I.  After  Fulton  (1977)  with  extrapolation  to 


J  osephson 
Junction 

Critical 
Current  Maximum 


Symbol 


Phase  Difference 


Voltages  x  2n/$ 


d<J)/dt 


polymer  dynamics 


Mechanical 


Pendulum 


Phase  Angle 
Pendulum 

Angular 

Velocity 


jCH^ 


u:  charge  dq  dt  related 
to  maximum  pseudo 
spin  vector  by 

dfl  =  U - 

dt  d<J>/dt 

=  gj  y  e J  /  dcf>  /  d  t 
d:  charge  dq^/dt 

u:  n— phase  kinking, 
d:  tt— phase  kinking. 

u:  voltage  x  2n/<J) 
related  to  spin 
angular  momentum  <o  = 
2pm  B/h  by  d$/dt  = 
hto/  (2m$  dq  /dt) . 
d:  voltage°x  5jt/$> 


•  o  .  • .  *  . '  .  •  «  -  ,-.-V  .  -V  • . ■ 


1 
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TABLE  I  (Cont'd) 


Josephson 

Janet ion  Symbol 

Applied  Current  I 


$o/(2nR)  D 


C$  / ( 2nR)  M 

o 


$  / (2nL)  K 

o 


Current  in  1^ 

Inductor 


Mechanical 
Applied  Torque 


Damping 

Coefficient 


Moment  of 
Inertia 


Torsion  Bar 
Spring  Constant 


Torsion  in 
Torsion  Bar 


!CHIx 

u:  applied  current, 
dq/dt,  related  to 
the  processional 
radius  of  the  spin 
vector  by  dq/dt  = 
#oj/(2m$  d$/dt) 

d:  appl ied  current , 
dq/dt. 

u:  none  for  ideal  sys¬ 
tem. 

d:  dq,/dx  =  D,  dissi¬ 
pation  to  external 
phonon  matrix. 

u:  da  /dx,  charge  re¬ 
pulsion. 

d:  da  /dx,  charge  re¬ 
pulsion. 

u:  dq^/dx,  excitonic 
coupling. 

d:  dqK/dx,  excitonic 
coupling. 

u:  dq^/dt  determining 
phonon  energy. 

d:  dq^/dt  determining 
phonon  energy. 


Josephson  Coupling 

Gravitational 

u 

Energy 

Potential  Energy 

-(  $  I  / (2n) )cos$ 

-I  cosd> 

o  c 

c  Y 

d 

Electrostatic 

Kinetic  Energy 

u 

Energy  1/2  CV* 

1/2  M(d<J>/dt)» 

-dqc/dt  cosc(>,  exci- 
ton-phonon  coupling 
energy. 

-dqc/dt  cos<|),  exci- 
ton-phonon  coupling 
energy. 

kinetic  energy  re¬ 
lated  to  paramagnet¬ 
ic  susceptibility  by 

(M)  *  _ 

Vdt ' 


1/2!5m 

dx 


X 

d: 


Jo 

2n 


d$ 

dt 

kinetic  energy: 
l/2dq^  ,dix  ’ 

ar  W 
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TABLE  I  (Cont'd) 


Josephson 

Inaction  Symbol 

Magnetic  Energy 
1/2  LI* 


Power  Input  from 
Applied  Current  IV 


Mechanical 


Torsional  Energy 
1/2I*/K 


Power  Input  from 
Applied  Torque, 
dq/dt.  d<f>/dt 


ICHix 

u:  1/2  (dq^) */dq^,  magn. 

dt  dx  energy 

d:  1/2  (dq^) */dq^,  magn. 

dt  dx  energy 

u:  dq/dt.  d$/dt,  power 
input  from  applied 
field. 

d:  dq/dt.  d$/dt,  power 
input  from  applied 
field. 


u  =  un doped  poly acetylene; 

d  =  doped  polyacetylene; 

«  =  h/2e  »  2.068  x  10-**  Webers; 

o 

and  Kd*<}>/dx*  *  I  sin<J>  +  Dd<j)/dt  +  Md*<|>/dt* 
c 


X 


dqT  44  =  dq  sin<()  dqn 

dt1*  dt  dt?dx  +  di° 

and  d$dx  =  LI.  *  ($  /2jr)d4>dx. 

Ls  O 


djfe 

dt 


+ 


+  i 


& 

dt* 


x’ 

+ 


or 


4a _ 

dt  dx  ' 


Some  predictions  can  be  made  for  the  current-voltage  relations  of  poly¬ 
acetylene  on  the  basis  of  this  model.  Firstly,  in  steady  state  there  is  an 
average  current  balance  given  by: 

<I>  =  <1  sin<J»  +  <V>/R.  (6) 

c 

Decomposing  I  and  V  into  steady  and  alternating  components:  I  =  <I>  +  <1  > 

and  V  =  <V>  +  <V  >,  then  if  I  is  fixed  <1  sin<J>>  =  <V*  >/R  and  there  should 
ac  c  ac 

be  a  departure  of  the  I-V  curve  from  ohmic  behavior,  indicating  that  the 

voltage  has  an  alternating  component. 

Secondly,  if  a  radio  frequency  current  at  frequency,  f,  is  applied  near 

the  junction  frequency  <V>/$  ,  the  ac  supercurrent  locks  into  phase 

o 

coherence  with  applied  currents  so  that  <V>  =  f  giving  linearity. 
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In  order  to  obtain  the  characteristic  particle  in  a  doable  potential 
($)  well  with  harmonic  spring  connection  to  neighbor(s)  (Figure  2), 

$ 

<t> 

Fig.  2 

an  inductively  connected  element  is  required  with  phases: $  =  2n/(V/$o>dt  and  @ 

■  2 n/(V  / $  )dt,  where  V  and  V  are  the  voltages  across  the  junction  and  an 
z  o  z 

external  passive  element  z,  i.e.,  the  next  neighbor  (Fulton,  1977): 


Fig.  3 


Alternatively,  the  <j>  of  Figure  2  is  the  angle  of  the  carbon  bond  defining 
the  isomerized  state  in  (CH)^,  and  the  $  is  the  configurational  potential. 
Thus,  the  two  potential  well  defines  two  isomerized  states  in  polymers. 

If  z,  the  external  passive  element  is  considered  a  resistor,  r,  (z  =  r), 
then  for  an  undamped  junction  having  LIc/$  >>  l/2n,  there  is  the  possibility 
of  a  relaxation  oscillation  if  L/ r  >>  (LC)1/*.  If,  on  the  other  hand,  z  is 
the  capacitor,  C  ,  and  we  define  C  .  =  (1/C  +  1/C  )-*,  then  a  resonance 

Z  6  X I  Z 

occurs  at  a  frequency:  co  =  (LC  __)1/a,  when  the  impedance  of  the  C-L-C 
combination  becomes  infinite. 

An  inquiry  may  now  be  made  concerning  what  happens  when  a  polymer  such 
as  polyacetylene  undergoes  an  order-disorder-order  transition  with 
displacement,  or  isomerized  state  transition,  e.g.,  for  one  element 
(monomer) : 


-10- 
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Fig.  4 

For  this  situation  of  the  monomer,  z  is  an  inductance.  L,  giving  a  T-bar 

twist.  LIq/$o,  balancing  the  pendulum  in  a  horizontal  position.  For 

LI  /$  <<  l/2n<  <J>  follows  6.  but  at  LI  /$  =  l/2n  the  pendulum  is  no  longer 

co  c 

stable  and  for  LI  /$>  l/2n  a  flip  to  the  other  state  occurs  (Fulton.  1977): 
c 


<£=3t/2 


Extrapolating  this  element  to  a  one-dimensional  chain  analogous  to  the 
continuous  limit  of  a  dc  SQUID,  the  equivalent  circuit  equation  is: 

K^d’cji/dx*  -  I  £  sin<(>  +  D£d4>/dt  +  M^d^/dt*  +  1^.  (7) 

where  IQ^  is  the  critical  current/unit  length,  1^  is  the  current/unit 
length  applied  to  the  junction,  *  $oI/2jr  q(2  +d)  =  L^  is  the 

inductance /unit  length,  =*  'J'0C^/2jt,  is  the  capacitance/ 

I  *  unit  length,  *  aW$o/2n,  and  aW  is  the  conductance/unit  length.  We  are 

especially  interested  in  the  behavior  of  this  chain  for  l  >>  Aj,  where  Aj  = 

I 

|  (<I,o/(2ffPo(2A+d)  jc> )  x/*,  when  vortex  propagation  is  exhibited.  Current 

!  propagation  is  described  by: 

! 

I 
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KjdH/dx*  -  Ic£sin4>  +  1^,  (8) 

which  is  in  the  same  form  as: 

Md^/dt*  -  -I  sin*}*1  +  I  ,  (9) 

c  x 

describing  the  motion  of  an  undamped  pendulum  subject  to  an  external  torque 
(see  equation  (5),  above).  A  vortex  (or  propagating  polymer  soli ton)  in 
this  analogy  is  as  follows:  a  pendulum  starts  in  a  pendulum  bob-down 
position  (order)  swings  up  (disorder)  falls  over  through  2n  and  comes  to 
rest  in  a  pendulum  bob-down  position  again  (order).  Conservation  of  kinetic 
and  potential  energy  gives  l/2(dip/dx)a  -  (1-co s<jj)AJ  for  such  motion, 
resulting  in  the  description  given  by  Figure  6.  The  two  bob-down  positions 
correspond  to  the  two  polymer  isomer ized  states. 

It  is  known  that  in  a  uniform  pendulum  array  of  infinite  length,  a  full 
twist  or  vortex  can  be  moved  from  place  to  place  with  no  cost  in  energy 
(Scott,  1969:  Waldram  et  al,  1970:  Lebwhol  and  Stephen,  1967).  Application 
of  a  weak  torque  in  the  vicinity  of  a  vortex  causes  the  pendulums  to  rotate 
in  the  direction  of  the  applied  torque  and  the  twist  moves  along  the  array 
away  from  the  torque.  This  is  the  propagating  soliton  of  the  displacive 
mode  of  the  system.  If  two  solitons  are  propagated  along  a  pendulum  array 
then  in  an  overdamped  system  the  two  will  be  annihilated,  in  an  underdamped 
system  the  two  pass  right  through  each  other  and  continue  on  with  their 
relative  kinetic  energies  intact. 

Consideration  may  now  be  given  concerning  how  it  is  that  the  vortex,  or 
n-phase-kink,  becomes  a  propagating  soliton.  Stated  differently,  how  does 
the  order-disorder  transition  mode  become  a  displacement  transition  mode? 

At  the  heart  of  the  matter  is  the  demonstration  by  Pitaevskii  (1959),  that 


for  a  system: 
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mm 


rig.  6.  Various  views  (a)  to  (d)  of  a  segment  of  an 
infinitely  long  pendulum  array  containing  a 
vortex,  (a):  Front  view  (cos<|)(x)),  (b):  Top 
view  (-sin  <J>  (x) ) ,  (c):  d4>/dx,  (d):  (x) .  (e): 

A  vortex  trapped  in  an  array  of  finite  length  by  the  excessive 
weight  of  the  end  pendulums.  (From  Fulton.  1977.) 
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the  change  of  the  complex  function  with  time,  t,  is  described  by  the 
equation:  ihBH'/dt  *  L  where  L  is  a  nonlinear  operator  and  a  function  of  a 

second  variable  —  temperature  in  the  system  considered  by  Pitaevskii.  As 
this  second  variable  increases,  the  anti-Hermitian  part  of  L  becomes 
dominant  (as  compared  with  the  Hermitian).  In  this  case,  the  equation  which 
describes  the  vortex  filament  goes  from  Schrodinger  type  to  diffusive  type 
(Pitaevskii,  1959,  1961:  Usui,  1969).  In  effect,  the  anti-Hermitian  part  of 
L  is  larger,  the  faster  the  system  relaxes.  L  cannot  be  Hermitian  for  a 
propagating  wave  (soliton),  because  density  would  then  be  conserved  and 

A 

relation  is  required.  If  L  is  Hermitian,  then  the  wave  is  static  (pinned). 

The  anti-Hermitian  part  describes  the  process  whereby  thermodynamic 
equilibrium  is  approached,  or  describes  the  process  whereby  the  density  of 
states  snd  monomer  1  in  a  chain  are  converted  into  the  density  of  states  at 
the  next  monomer,  monomer  2  (Barrett,  1983).  A  nonlinear  operator  gives: 

3E/3t  +  divQ  »  0,  wave  Q  is  an  energy  loss  vector.  For 
* ' superconductivity ' '  or  decrementless  conduction  between  monomers  1  and  2: 

divQx  +  divdj  -  0;  (10A) 

3E1/3t  +  3E2/3t  -  0.  (10B) 

A 

In  the  case  of  the  anharmonic  oscillator,  we  have  3<J>/3t  =  LV,  or 

<)>  »  C  n^n  (I  (X1  +  ix  +  C)  +  (u(^)  -  u(f/))»)  (IOC) 

2  U'  l _ B _  _ 

Hermitian  Anti-Hermitian 

Component  Component 

where  x  is  the  position  of,  e.g.,  a  carbon  bond  (for  (CH^)), 

<j>  is  configurational  potential,  and 
l  is  the  length  of  the  polymer. 

In  the  case  of  the  Davydov  vibrational  (protein)  soliton,  we  have 
3a/ 3t  »  La,  or 
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(i)  iXdft 
dt 


LULL 

Hermitian 


Component 

+  x,(r,  -  p  ,  )  a  -  J(a 

1  *^n+l,ti  *n-l,r 


no 


,  +  a  ,  )  + 

n-l,a  n+l,a 


+  L  Un.a+1  +  *n,a-l>  + 


+  *2^n+l,a  an+l,a  ^n-l,a  *n-l,a  ^na^*n+l,a  *n-l,a 


)) 


Anti-Hermitian  Component 


-  *(S0+1.a  -  ^„a  +  K-I.a' 
_&£_ - - - 


(10E) 


Hermitian  Component 

-  Vl‘.«l‘  -  IVil‘> 

Anti-Hermitian  Component 


where  and  X^  are  nonlinear  coupling  coefficients. 


T  and  L  are  dipole  coupling  coefficients. 


E  ■  1650  car-1  ■  0.205  e.v., 
o 


a  is  the  probability  of  finding  a  quantum  of  bond  energy, 

8  is  longitudinal  displacement. 

Thus,  in  considering  whether  a  polymer  can  exhibit  soliton  behavior, 
neither  the  existence  of  isomerized  states  nor  the  existence  of  vibrational 
modes  of  energy  compatible  with  an  energy  input  is  sufficient  condition  for 
propagating  waves.  Only  if  the  anti— Hermitian  component  of  the  equation  of 
state  exceeds  the  Hermitian  component  will  a  perturbation  become  dispersive 
or  propagate  (Barrett,  1983). 

Given  that  a  polymer  chain  will  support  a  soliton,  consideration  may  be 
given  as  to  whether  a  bifurcating  chain  (T-junction)  will  allow  soliton 
propagation  down  the  arms  of  the  chain. 


III.  T-JUNCTIONS 

Consider,  then,  three  chain  polymers  joined  at  one  end  to  form  a  T  or  I 
junction.  Such  joined  chains  may  be  referred  to  as  T  junctions.  In  the 
Josephson  line  of  one  dimension  Equ.  (2),  above,  becomes  (Scott,  1970): 
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1  0ft  *  0 ' 0  '  r  ft  "  *i0*  ‘  Y  (11) 

where  I  -  L(2nJ  /$C)VVr. 
c 

r  -  g($  lira  c>*/», 

o  c 

*  -vj. 

$  -  h/2e 

o 

J  m  a  constant  giving  the  maximum  Josephson  current  per  unit  length: 
c 

L,  C,  r,  g  and  jg  are  the  series  inductance,  the  shunt 

capacitance,  the  series  resistance,  the  shunt  conductance,  the 

distributed  bias,  and  the  current  source  per  unit  length,  and 

where  distances  are  measured  in  units  of  *  ($  C/2nJ  )*/*. 

J  o 

Simply  stated  K  and  T  are  losses  in  the  line  and  Y  represents 
bias. 

In  the  case  of  pendnla  line  T-junctions,  Naka j ima  et  al.  (1976)  define  two 
classes:  (1)  the  TTP  (trigger  turning  point),  which  is  analyzed  by  the  two- 

dimensional  equation  reported  by  Naka j ima  et  al.  (1974a  and  1978c).  The 
boundary  condition  for  the  TTP  is: 

<f>  |  ■  <f>  |  ■  <J>  |  (12) 

*1 '  tp  *2  tp  *3  *  tp 

(2)  A  second  turning  point  is  the  STP  (selective  turning  point),  with: 

4> 

xi  +  x2  +  x3  -  0.  (13) 

Thus,  a  single  flux  quantum  propagating  toward  an  STP  on  one  line  will 
initiate  a  single  flux  quantum  on  only  one  connected  line,  depending  on  the 
bias  current  of  each  line,  the  applied  magnetic  field  and  the  junction 
geometry. 

Confirming  attention  to  the  TTP  junction  and  commencing  with  the  two- 
dimensional  description  of  Equ.  (10),  above:  the  motion  of  a  vortex  filament 
at  temperatures  close  to  the  ^-point  may  be  described.  As  the  stability  of 
a  vortex  is  not  guaranteed  for  a  real  field  in  two  dimensions,  Nakajima  et 


al.  (1978c)  introduce  a  complex  field: 
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\p  *  Re  ^  + 

ilmY 

(14) 

with  now: 

(^ 

+  ^ 

ft  >  *"*'  *  9R# 

(15) 

3xa 

3y* 

(li_ 

l3x» 

+  ai- 

3y» 

- 

It  >  ^  *  IS* 

(16) 

where  the 

potentials 

UR 

and  Uj  are  related  by: 

V 

}  0M>)4  - 

1 

2 

[l-(Im^)*](R#)», 

(17) 

UI- 

J  (Imi{>)«  ■ 

1 

'  2 

[l-(Re^)*](Im^)». 

As  the  complex  field  possess  a  phase  angle,  clockwise  and  anti-clockwise 
vortices  are  obtained.  With  the  condition  stated  by  Equ.  (12)  satisfied, 
then,  nenristor-like  activities  are  possible  (Crane  (1962))  and  describable 
in  two-dimensional  form  (Eqn.  (10))  with  a  complex  field. 

The  problem  of  defining  specifications  for  propagation  of  solitons 

through  (a)  a  TIP  junction  and  (b)  STP  junction  may  be  approached  as 

follows:  In  the  case  of  (a),  the  TTP  junction,  the  requirement  is  given  in 

Equ.  (12).  As  <p(t)  *  (2n$  )jV( u)du,  we  require  the  voltage  x  2n/$  to  be 

o  *  o  o 

the  same  on  all  3  arms  of  the  T  junction.  This  reduces  to  the  simple 
conclusion  that  all  three  arms  must  be  of  similar  polymer  material  and  of 
the  same  length. 

Case  (b),  the  STP  junction,  is  more  complex.  Referring  now  to  Equ. 
(11),  for  a  set  T,  a  soliton  at  a  turning  point  may  (1)  halt,  (2)  propagate 
into  arm  1,  (3)  propagate  into  arm  2,  or  (4)  propagate  into  arm  3,  depending 
upon  the  ratio  of  the  bias  currents:  Y  corresP°nds  to  the  applied 

torque,  I,  of  the  single  junction  (Table  I).  This,  in  turn,  is  a  function 
of  the  amount  of  "pinning",  i.e.,  introduction  of  defects  into  the 
polyacetylene,  in  order  to  obtain  critical  current  increases  at  an 
appropriate  applied  field  (Jaklevic  et  al.  1965,  Tamashita  et  al.  1974).  In 
order  to  predict  which  of  the  behaviors  (1)  -  (4)  will  occur,  one  also  needs 
an  estimation  of  d/2Aj,  where  2^j  is  the  length  of  the  soliton  and  d  is  the 
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length  of  the  polyacetylene  (Owen  and  Scalapino,  1967).  An  easy  conceptual 
solution  is  thus  to  uniformly  dope  a  length  of  (CH)^  (i.e.,  set  Y^  =  = 

Yj)  and  synthesize  T-junctions  of  different  arm  lengths  for  different 


Fig.  7 

In  the  vicinity  of  the  junction  J,  the  normal  phonon  modes  will  be 
different  from  those  at  the  ends  Mid.  £3*  o*  in  the  center  of  each 

chain.  As  soliton  propagation  of  the  moving  charge  domain-wall-type  is 
dependent  upon  exc it on-phonon  coupling,  the  arrival  of  a  soliton  in  the 
vicinity  of  J  should  show  different  propagating  behavior  compared  with  its 
behavior  at  a  distance  from  the  junction,  due  to  different  phonon-dependent 
propagation  energies  for  that  carbon  atom  at  the  T  junction  marked  with  a  J 
(Fig.  7).  This  situation  is  different,  therefore,  from  the  pendulum  line  as 
the  T  junction  in  polymeric  form  exhibits  different  monomer  vibrational 
forms  at  the  T  join;  and  the  likelihood  is  that  the  soliton  would  not 
progress  past  the  junction  into  arms  d^  and  d^  but  be  pinned  there. 

IV.  TUNNELING 

In  the  previous  sections,  treatment  of  soliton  propagation  in  a  one¬ 
dimensional  organic  conductor,  trans-polyacetylene,  was  made  in  analogy  with 
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fluxon  propagation.  This  treatment  is  extended  here  to  the  case  of  uneven 
doping  and  unevenly  conjugated  conductors. 

Parametric  excitation  is  the  underlying  dynamic  of  the  pendulum  analogy 
as  well  as  the  Josephson  effect. 

fe  commence  with  Eq.  (3)  of  section  II  for  the  total  free  energy  for  a 
Josephson  barrier: 


p  /dx[(l-cos<J>)  +  i  <1^  >a  +  |  <5  ) 

2e  2  o  ox  2  c 

in  which  the  rate  of  energy  dissipation  is: 

■  <1“co,<t,)  +  2  Xo(l^)a  +  2  (r-)Mlt)a* 


Recalling  that  3<j>/3t  is  the  mechanical  angular  velocity  or  the  Josephson 

junction  voltage  x  2n/<J>0#  we  now  extend  the  voltage  analogy  to  3<J>/3x  and 

call  this  term  the  dispersion  voltage.  As  (l-cos<j>)  is  the  coupling  energy, 

it  is  evident  that  Eq.  (19)  satisfies  the  condition  for  parametric 

excitation,  where  the  l^r11) * (“^) *  term  is  the  pump,  (l-cos<j>)  is  the  signal 

and  1/2  is  the  idler.  Regarding  a  Josephson  junction  or  polymer 

o  <3x 

chain  from  a  parametric  excitation  point  of  view  permits  the  treatment  of 
the  line  in  discrete  elements,  dx.  Variation  of  d<f>/dt  within  any  dx  results 
in  a  variation  of  d<j>/dx,  or  dispersive  transmission,  due  to  the  conservation 
condition  of  Eq.  (19).  The  monomer  for  the  polymer  chain  is  described  by 
the  equation: 

<p"  m  +  af  +  b  =  0,  (20) 

o  o 

where  fQ  siay  be  taken  as  the  midfrequency  of  a  signal  equivalent  to  electron 
momentum  in  the  analysis  of  states.  Eq.  (20)  is  formally  equivalent  to 
Schr'ddinger's  equation.  If  the  following  substitutions  are  made: 
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6uu 


(21) 


1X1 


‘V 


\  ■*  b' 


where  subscripts  indicate  partial  differentiation,  then  one  has  the  Korteweg 
de  Vries  equation  (Korteweg  and  de  Vries,  189S): 


ut  +  6uu^  + 


u  =0 
xxx 


which  is  derived  from  the  Boussinesq  equation  (Boussinesq,  1871): 


tt 


r-  (u  +  6uu  +  n  ), 
dx  z  z  zzz 


(22) 


(23) 


by  restricting  the  motion  to  one  direction.  Ve  note  that  Eqs.  (20)  and  (22) 
are  formally  similar. 

If  now,  we  consider  a  light  wave  polarizing  energy  levels  of  a  polymer 
with,  for  example,  two  levels  of  energy,  then  the  polarization  may  be 
decomposed  into  components  in  phase  and  in  quadrature  with  the  electromag¬ 
netic  wave.  The  Schrodinger-type  equation  (20)  for  the  atom  reduces  then  to 
the  Bloch  equations  for  polarization  or  the  self-induced  transparency 
equations  (Lamb,  1971).  After  transformation,  these  equations  yield  (Scott 
et  al.,  1973): 


«u  “  utt  *  sinu,  (24) 

or  the  sine-Gordon  equation. 

The  sine-Gordon  equation  is  the  unique  equation  invariant  under  the 
infinitesimal  Lie  transformation  z-*"  (1  -e  )x,  t"*"  (1  +  £)t  and  the  finite 
scale  transformation  a**  a-1z,  f*  at  and  is  Lorentz  invariant  (Bui lough  and 
Caudrey,  1980) . 

Now,  the  relations  governing  power  flow  in  lossless  nonlinear  lumped 
circuit  elements  are  the  Manley-Rowe  relations  (Manley  and  Rowe,  1956). 
Considering  the  mean  flux,  N^,  of  photons  of  frequency  u,  the  analogous 
relations  for  extended  lossless  nonlinear  media  are  (McLean,  1977): 
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Nul  +  Nu2  +  Nu  3  =  constant' 


Q1  +  u2  =  °3  * 


Substituting: 


u  ->  u, , 

XX  1 


Utt  V 


sm  u  -*■  u^. 


provides  a  formal  similarity  between  the  sine-Gordon  eq.  (24)  and  the 
equation  describing  parametric  amplification  (Eq.  (26)).  The  relation  of 
the  sine— Gordon  equation  to  electron  tunnelling  and  the  Josephson  effect  is 
already  known  (Bullough,  1977),  by  which  the  parallelism  is  drawn  between 
the  ''kink''  and  ''anti-kink''  solutions  of  the  sine-Gordon  equation  and  the 
fluxons  in  large  area  Josephson  junctions.  Here,  we  point  out  the  parallels 
between:  (1)  the  Schrodinger  formalism  and  the  Eorteweg-de  Vries  equation 
and  (2)  the  sine-Gordon  equation  and  the  Manley-Rowe  conservation  equations. 

We  now  utilize  this  development.  Suppose  the  wave  function  for  an 
electron,  i1,  is  given  as  a  function  of  distance,  x,  the  Schrodinger  equation 
for  which  is: 

-  &  ^  *  v<‘>*  ■  E  •  <28> 

where  V(x)  is  the  amplitude  of  a  barrier  to  free  movement  of  the  electron. 

If  V  varies  relatively  slowly  with  distance,  then  the  solution  of  Eq.  (28) 
is  (Messiah,  1958): 

iJj  =  ip  exp[±i  Is  /E  -  V(x)  x]  (29) 

o  h 

‘  *o  eiI,ttl(/<"p»p  -  uSig0.l,),!- 
Adapting  Eq.  (26)  to  explicit  parametric  amplification  form,  we  have: 


u  =  u  .  ,  +  u.., 

pump  signal  idler 
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Eqs.  (29)  and  (30)  are  then  interpreted  as  follows:  ^7  E  and 

^7  V(x)  are  the  spatial  frequencies  of  two  oscillating  systems,  Eq.  (29) 

represents  the  situation  of  parametric  excitation  in  which  7®  E  =  u  is 

h*  pump 

the  spatial  frequency  of  the  pump,  ^7.  V(x)  =  u  .  ,  is  the  spatial 

ft*  signal 

2m 

frequency  of  the  signal,  and  ^7  (E-V(x))  =  u^j  is  the  spatial  frequency 
of  the  idler.  As  (i),  by  definition,  V(x)  is  maximum  at  a  different 
position,  x,  than  is  i.e.,  it  is  a  barrier,  and  as  (ii)  it  is  the  pumped 
signal  system  which  is  of  energy  V(x),  it  follows  that  a  pumped  idler  system 
of  spatial  frequency  (E-V(x))  should  appear  at  a  distance  x,  i.e.,  out  of 
the  barrier.  Vhereas  Joseph son  (1962)  treated  the  term  in  the  Hamiltonian 
which  transfers  electrons  across  the  barrier  as  a  perturbation,  here  we 
specifically  designate  it  as  a  pumping  system.  Electron  tunnelling  is  thus 
displacive  parametric  excitation  and  may  be  given  a  classical  mechanics 
point  of  view. 

For  example,  the  master  equation  for  the  effect  is: 

<31> 

In  the  case  of  two  superconductors  connected  by  a  ''weak  link'',  the 
Josephson  equations  (Josephson,  1962,  1964,  1974)  give: 
j  -  -  *!  -  {S  dx). 

It  (*2  -  *1  - 1  d”  -  -  i  <ui  -  u2>  -  f*  '  ‘«> 

for  a  one-dimensional  junction,  and  where  j  is  the  supercurrent  through  the 
junction,  4>  is  the  phase  of  the  order  parameter,  A  is  the  magnetic  vector 
potential,  U  is  the  chemical  potential,  and  V  is  the  voltage  bias  across  the 
junction.  Again,  the  junction  bias  acts  as  a  signal  to  the  second  system  as 
pump  and  the  third  as  idler. 


This  is  apparent  in  Josephson's  (1964)  account  of  the  tunnelling 
effect.  As  the  barrier  is  changed  from  a  very  thick  one  to  a  vanishingly 
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thin  one,  the  system  loses  the  degree  of  freedom  corresponding  to  the 
ability  to  alter  phases  independently.  The  free  energy  of  the  system 
contains  a  contribntion  from  the  barrier  region  which  depends  on  the 
relative  phases  of  the  values  of  $  on  the  two  sides  of  the  barrier,  and 
whose  magnitude  becomes  greater  as  the  barrier  is  made  thinner.  Perfect 
conductivity  of  a  superconductor  is  the  consequence  of  the  fact  that  in  a 
steady-state  the  order  parameter  must  have  the  form  (Josephson,  1964): 

y(r.  t)  =  VU,  0)  exp  [i  a(t)3  .  (33) 

I'  will  still  have  this  form  in  a  system  partitioned  by  a  barrier,  provided 
that  the  phases  of  the  order  parameters  on  each  side  of  the  barrier  are 
locked  together  by  the  barrier  energy.  The  only  way  this  can  occur  is  by 
parametric  excitation. 

V.  EFFECT  OF  PINNING 

In  sections  I  -  III  mention  was  made  of  the  effect  of  removal  of 

pendula  for  a  length  larger  than  2\  ,  where  A  is  the  Josephson  penetration 

J  J 

length  or  O^/mgl)1/*  in  the  model  (Kq  is  the  twisting  constant  of  the 
elastic  coupling  per  unit  length,  1  is  the  pendulum  length,  mg  is 
gravitation),  resulting  in  the  pinning  effect  of  a  Josephson  junction. 
Referring  to  Eq.  (19)  above,  the  effect  of  pinning  results  in  the  neglect  of 
the  (l-cos6)  term,  i.e.,  0-K).  This  term  is  the  signal  term  of  the 
parametric  excitation  dynamic,  hence  there  can  be  no  conduction  beyond  such 
a  barrier.  Conduction  between  barriers  along  a  line  of  intact  pendula  is, 
however,  possible.  Such  a  segment  would  exhibit  the  properties  of  a  soliton 
if  short  time  constants  are  considered,  and  a  trapped  polaron  if  a  long  time 
constant  is  considered. 

Pinning  may  be  obtained  in  another  way.  The  Josephson  penetration 
length  is  a  function  of  mg,  hence  weighting  the  pendula  bobs  would  obtain 
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the  same  effect  as  removal  of  pendula.  A  pinned  soliton  (kink)  is  shown 
below: 


Fig.  8 

The  result  of  pinning  (or  doping)  is  firstly  to  provide  a  barrier  to 
conduction.  However,  the  critical  Josephson  current  increases  at  an 
appropriate  applied  field.  From  the  mechanical  point  of  view  this  means 
that  a  torque  is  placed  on  the  first  pendulum  of  a  magnitude  such  that 
either  (a)  the  barrier  of  the  vacated  pendulum  is  traversed,  or  (b)  the 
heavy  bob  of  the  second  pendulum  is  raised.  Either  way,  an  angle  ®  is 
created  for  the  second  pendulum  permitting  a  signal  term  (l-cos6)  for 
parametric  excitation. 

The  effect  of  nonuniform  doping  in,  for  example,  (CH)^,  mav  be 
considered,  i.e.,  when  pendula  are  removed  so  that  irregular  spacings  are 
created,  or  when  the  pendula  have  bobs  of  varying  weight.  These  variations 
all  result  in  variations  in  the  9s  created  at  each  dx  along  the  chain. 

This,  in  turn,  results  in  variations  in  the  (l-cos9)  or  signal  term  at  each 
dx.  The  conservation  condition  of  parametric  excitation  then  requires  that 
increased  doping,  resulting  in  a  smaller  6  and  hence  a  smaller  signal  term, 
gives  a  larger  idler  voltage.  In  short,  increased  doping  should  result  in  a 
larger  voltage  along  the  line  due  to  increased  resistance  at  constant 
current.  Turned  around,  at  constant  voltage  increased  doping  results  in 
increased  resistance  hence  decreased  current.  Simply  stated,  the  tolerance 
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of  as  irregularly  doped  trans-polyacetylene  chain  is  that  of  a  line  with 
resistances  of  various  values  in  series. 

Considering  the  problem  from  a  more  basic  point  of  view.  Eq.  (20) 


above,  is  of  the  same  form  as: 


3*¥  3 *¥ 

3x*  +  3y* 


3T 

3t 


(|V|*  -  1)'}', 


(34) 


which  is  Eq.  (10)  of  section  II.  There,  the  demonstration  by  Pitaevskii 
(1959)  was  mentioned  that  for  a  system  such  as  this,  the  change  of  a  complex 

A  A 

function  I*  with  time  is  described  by  the  equation:  ihS'F/St^L'y,  where  L  is  a 
nonlinear  operator.  As  the  second  variable  increases,  the  anti-Hermitian 
part  of  L  becomes  dominant  (as  compared  with  the  Hermitian) .  In  such  an 
event,  the  equation  which  describes  the  polymer  chain  goes  from  Schrodinger 
type  to  diffusive  type  (Pitaevskii,  1959,  1961;  Usui,  1969).  Thus,  the 
effect  of  an  increase  in  doping,  beyond  a  certain  level  discussed  below,  is 
an  increase  in  Schrodinger  type  behavior,  the  effect  of  a  decrease  in 
doping,  below  this  level,  is  to  increase  again  diffusive  type  behavior.  For 
a  room  temperature  type  conductor  which  is  uniformly  and  ideally  doped,  we 
could  obtain  the  same  effect  over  time  of  nonuniform  doping  by  a  temperature 
fluctuation  which  resulted  in  a  variation  in  the  ratio  of  Schrodinger  and 
diffusive  type  behavior.  When  there  is  behavior  described  only  by  a 
Schrodinger  type  equation,  there  may  be  a  trapped  polaron,  but  there  is  no 
conduction  beyond  the  length  considered. 

VI.  OHM'S  AND  gIRCHOFF ' S  LAWS 

In  recapitulation:  in  a  polymer  chain  of  monomer  length,  dx,  there  can 
be  either  order-disorder  Schrodinger  type  behavior  or  displacive  behavior. 

In  an  electrically  conducting  polymer,  the  former  behavior  results  in 
resistance,  the  latter,  in  conductance.  These  two  modes  can  be  related  to 
the  conservation  laws:  Olm's  law  and  Kirchoff's  law. 
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Halperin,  Hohenberg  and  Ma  (1974,  1976)  aad  Hohenberg  and  Halperin 
(1977)  have  studied  the  effects  of  different  conservations  on  the  energy  and 
space  integral  of  the  order  parameter.  The  order  parameter  in  (CH)^  is  4>  or 
the  pseudo  spin  (undoped  case)  or  bound  charge  (doped  case),  which  can 
exchange  energy  with  a  reservoir  that  has  either  an  infinite  specific  heat 
or  an  infinite  thermal  conductivity.  The  first  case  typifies  resistance 
conduction  and  the  second,  decrementless  conduction.  Energy  conservation 
can  be  achieved  by  choosing  the  damping  constant  r  appropriately.  Four 
cases  are  then  distinguished:  (a)  <t>  not  conserved,  E  not  conserved,  and 

is  a  constant;  (b)  <j>  conserved,  E  not  conserved,  T  *  -  X  V;  (c)  $  not 

o  o 

conserved,  E  conserved, T  is  a  constant;  (d)  <p  conserved,  E  conserved, 

o 

r  -  -  XQV.  In  terms  of  the  mechanical  analog,  conservation  of  the  pseudo 

spin  corresponds  to  conservation  of4>:  conservation  of  energy  corresponds  to 

conservation  of  the  applied  torque;  the  constant  To  corresponds  to  the 

damping  coefficient  D  and  moment  of  inertia  N  of  Table  I  in  section  I;  thus, 

T  *  -  X  V  indicates  a  known  variation  in  the  size  and  weight  of  pendula 
o  o 

bobs. 

In  the  analogy  drawn,  conservation  of  <J>  will  occur  only  if  Kirchoff's 
law  applies  to  the  system  as  an  adiabatic  system.  Furthermore,  Obn's  law 
entails  conservation  of  energy  in  the  system  considered  adiabatic.  The  four 
cases,  A-D  may  thus  be  represented  in  terms  of  four  nonadiabatic  circuits 
(heavy  lines)  joined  to  reservoir  circuits  in  series  (lighter  lines).  The 
elements  of  the  circuits  consist  of  a  battery,  resistances  and  lead  to 
ground: 
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Mechanical 


Trans  (CH) 


T  junction  1. 

conpling  to 
other  pendulum 
line. 

T  junction  2. 

coupling  to 
other  elastic 
torque  line. 

Pendula  bob  3. 

weights  on 
original  line  equal. 


T  junction  coupling 
to  other  fibers 
unequally  doped. 

T  junction  coupling 
to  other  fibers  of 
unequal  ezcitonic 
coupling. 

Original  fiber 
uniformly  doped. 


No  T  junction  1. 

to  other  pendulum 
line. 

T  junction  to  2, 

line  offering 
elastic  torque 
only. 

Pendula  bob  3. 

weights  on 
original  line 
unequal. 


No  T  junction 


T  junction  to  ligand 
offering  different 
ezcitonic  inter¬ 
action  from  original 
fiber . 

Original  fiber  non¬ 
un  if ormly  doped. 


Mechanical 


Trans  (CH) 


T  junction  1. 

coupling  to  other 
pendulum  line. 

No  T  junction  to  2. 
line  offering 
elastic  torque. 
Pendula  bob  3. 

weights  on 
original  line  equal. 


T  junction  coupling 
to  fibers  unequally 
doped. 

No  T  junction  to 
fibers  offering 
ezcitonic  coupling. 
Original  fiber 
uniformly  doped. 


No  T  junction 
to  other  pendulum 
line . 

No  T  junction 
to  other  pendulum 
1  ine . 

Pendulum  bob 
weights  on 
original  line 
unequal. 


1.  No  T  junction 


2.  No  T  junction 


Original  fiber 
nonun  if ormly  doped 
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This  description  of  trans-polyacetylene  only  applies  to  the  displacive 
node.  This  nay  be  defined  more  accurately  by  consideration  of  the  following 
classical  Hamiltonian: 


H 


I  (•£*— 

i  2Mo 


*2  Vi 


W**- V*' 


(35) 


where  the  subscripts  i  and  j  refer  to  points  on  a  d-dimensional,  simple 
cubic  lattice  with  lattice  constant  unity,  T  is  a  scalar  quantity 
describing  the  displacement  of  an  atom  in  the  ith  unit  cell,  and  p^  = 
t^dH^/dt  is  the  momentum  conjugate  to  The  first  term  is  of  the  same 

form  as  Eqs.  (20)  and  (34)  above.  If  | rQ [ / Jc  (>>1  the  transition  undergone  by 
the  system  is  an  order-disorder  (Schrodinger-like)  transition:  if 
|ro(/Jo>>l,  the  system  undergoes  a  displacive  transition  (Hohenberg  and 
Halperin,  1977).  The  above  description  for  A,  B,  C,  and  D  refers  to  the 
displacive  transition  only. 

VII.  DYNAMIC  FORM  FACTOR 


Conducting  trans-polyacetylene,  proteins  and  DNA  with  soliton 
excitation  are  in  the  displacive  transition  mode,  which  exhibits  unique 
attributes.  The  critical  dynamics  of  such  a  phase  transition  is  dominated 
by:  (1)  a  soft  phonon  whose  frequency  decreases  as  the  transition 
temperature  is  approached,  and  (2)  the  appearance  of  a  central  component 
around  w=0  in  the  dynamic  form  factor  (Schneider  and  Stoll,  1973,  1976, 
1978,  1981).  This  dynamic  form  factor  or  spectral  density  is  related  to, 
and  defined  by,  the  characteristic  frequency  or  relaxation  time  of  a  mode: 


'A*  <k> 


t ..(ki  -  s  (k,^o> .  2 /“  fAAj , 

AA  AA  o  S  (k,t=0) 


where  the  correlation  functions  are  (Schneider  and  Stoll,  1976) 
1 


SAA(k,t)  B  SAA(k,t=0)  “AA 


SiA(k,t),  S^(k,t)  =  <A(-k,0)A(k,  t)  > , 


and  related  to  the  spectra  densities  or  dynamic  form  factor  by: 
2 

J  f*  ac  mr  ' 

AA 


SAA<1'*)  -  S“  (k-t,dt- 


(36) 


(37) 


(38) 


AA 
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At  the  present  time  there  is  no  firmly  based  analytical  technique  to 
calculate  the  quantities  (36) ,  (37)  and  (38)  for  a  system  evolving  according 
to  Newton's  equations  (Schneider  and  Stoll,  1981).  These  quantities  are, 
however,  of  crucial  importance  in  understanding  the  excitation  spectrum. 
Schneider  and  Stoll  (1973)  reported  that  during  the  temporal  development  of 
this  system  clusters  form,  representing  connected  cells  whose  instantaneous 
local  order  parameter  has  a  sign  opposite  to  that  expected.  These 
propagating  clusters  in  two  or  three  dimensions  are  the  equivalent  of 
soliton  propagation  in  the  one -dimensional  trans-polyacetylene.  For  small 
wave  vectors  where  the  frequency  of  the  cluster  waves  is  small,  the  damped 
cluster  waves  become  overdamped  and  give  rise  to  the  central  peak 
phenomenon.  The  central  peak  half-width  is  proportional  to  the  inverse 
lifetime  of  the  cluster  waves  and  its  height  is  proportional  to  this 
lifetime  (Schneider  and  Stoll,  1976).  The  behavior  of  these  clusters  is 
consistent  with  the  universality  hypothesis  for  structural  phase 
transitions.  The  prediction  may  be  made,  then,  that  (a)  neutron  scattering 
will  detect  a  central  component  at  w-0  in  the  dynamic  form  factor  in  the 
case  of  polymers  excited  to  soliton  propagation;  and  (b)  that  the  dynamics 
of  such  polymers  obey  the  universality  hypothesis,  if  soliton  propagation 
occurs. 

This  development  can  be  placed  in  the  context  of  the  thrust  of  sections 
I-III.  Iwamatsu  and  Onodera  (1977)  have  derived  an  exact  expression  for  the 
dynamical  structure  factor  within  the  framework  of  classical  statistical 
mechanics.  These  investigators  commence  with  an  equation  analogous  to  Eq. 
(20): 


H  =  M  x*  +  Ax4  +  Bxa 
2 
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where  M  is  the  ness  of  the  oscillator,  z  is  the  displacement  of  the 
oscillator  from  the  origin  or  normal  coordinate  of  a  particular  mode,  B 
represents  the  coefficient  of  harmonic  restoring  force  and  A  represents  the 
contribntion  of  enharmonic ity.  The  anharmonic  potential  7(x)  *  Ax4  +  Bx* 
has  one  or  two  minima  depending  on  the  sign  of  B: 

i  ■fcm 


Fig.  10 
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However,  the  results  obtained  can  be  generalized  to  a  nonlinear  oscillator 
moving  in  a  quite  general  potential  (Matsubara,  1972),  for  example: 


Fig.  11  \lgt 

An  equivalence  between  Eq.  (39)  and  Eq.  (20)  is  obtained  for:  a=B/A, 

b-Mx1/ 2A,  B=aA,  and  the  energy  E  is  in  units  of  B*/4A  *  a*A/4.  In  the  case 

of  E<0  (B<0)  the  motion  of  the  oscillator  is  confined  within  one  of  two 

wells  and  the  solution  is  the  Jacobian  elliptic  function  (Onodera,  1970): 

x(t)  *  adn(ft  t  -  dn-i(x/a))  (40] 

o 

where  ■  (1+E)1/4.  The  fundamental  frequency  of  oscillation  is: 


2jfto  m  _JL_ 
2K(k)  K(k) 


+  1  +  E)V» 
2 


where  K(k)  is  the  complete  elliptic  integral  of  the  first  kind.  The 
dependence  of  0(E)  on  the  energy  E  is  (Iwamatsu  and  Onodera,  1977): 


of-  Y'b.j 


0  i  2  J  « 


Fig.  12 
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The  degree  of  enharmonic ity  is  given  by  the  modulus  k: 

k*  -  1/2  (1  -  Tfijf),  B  >  0,  (■ 

k*  =  1  +  (1+E)-1/1  B  <  °»  (- 

and  the  amplitude  of  harmonic  oscillation  and  anharmonic  oscillation  is: 

a*  -  2(/l+E  -  sgnB) ,  B  >  0,  (. 

a*  -  2(1  +  /l+E),  B  <  0,  <- 

a1  =  J  (1  -  /lHc4) .  (, 

The  oscillation  amplitude  versus  energy  for  tiro  wells  is: 


5  -  -J  I  I  3  l  ! 


Fig.  13  ,w 


The  dynamic  structure  factor  S(k,w)  is: 

S(k,w)  *  S^(k,w)/Z,  for  single-minimum  potential  case  (44 

S(k«w)  =  ^2^*+^2Z^  ^kl for  the  double-minimum  potential  case  (44 

where  Z,  -  f°  dEe“E/T  /ME)  . 

1  o  ' 

S1(k.w)  -  /”  dEe_E/T^jB_00 1 a1  (k,E)  I  *  d(w)-nH(E)), 

n 

ot j  (k,E)  is  the  Fourier  coefficient  of  the  periodic  function  exp(ikacn^  t) 
n  0 


n  /V  FI  _  1  ,T(E)  ikacnft  t  -inft(E)t.. 

ai  (k*E)  *  T(E)  fo  9  °  e  dt' 

n 


S2(k.w)  -/"j  dE.  E/T  (k,E)  |  * 


Z,  -  /!2dE.-E/I  .»d 

a2  is  the  Fourier  coefficient  of  the  periodic  function  exp(ikadn^  t) 

n  ° 


a-  (k.E) 


1  rT(E)  ikadnH  t  -infi(E)t. 
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It  follows  from  these  expressions  that  la^nl  and  jc^  |  have  local  maxima 
when  the  transferred  wavelength  A=2n/k  or  its  positive  integer  multiples,  is 
nearly  equal  to  the  oscillation  amplitude,  i.e.,  when: 

ka  or  ka1  ~  a  (2p-l),  (p=l,2,...).  (45) 

2 

This  means  that  particular  oscillations  that  have  a  oscillation  amplitude 
comparable  with  the  transferred  (transmitted)  wavelength  (or  positive 
multiples  thereof)  may  resonate  with  a  perturbation  offered  by  a  neutron 
beam,  i.e.,  inelastic  neutron  scattering  is  a  possible  investigative  tool. 
This  has,  indeed,  been  demonstrated  for  SrTiO^,  EMnF^  and  KTaO^  (Shirane, 
1974).  The  expectation  is  that  trans-polyacetylene  will  demonstrate  similar 
scattering. 

The  dynamic  structural  phase  transition  accompanied  by  the  softening 
(i.e.,  by  acting  as  pump)  of  a  particular  mode  of  lattice  vibrations  is 
shown  in  the  following  example  (Iwamatsu  and  Onodera,  1977): 


0  I  Z  3  4 


Fig.  14 

This  constitutes  the  strongest  piece  of  evidence  that  soliton  propagation 
for  the  (CH)^  type  of  polymers  is  pumped  by  a  phonon  mode  with  phonon- 
exciton  coupling  as  stated  above. 

The  question  of  tolerance  in  soliton  propagation  in  an  irregularly 
doped  transpolyacetylene  line  thus  may  be  translated  into  a  calculation  of 
the  dynamic  form  factor  for  a  particular  nonlinear  oscillator  moving  in  a 
general  potential  at  a  particular  dx  on  the  line. 
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One  may  ask  what  is  so  special  about  organic  conductors?  Other  systems 
undergo  phase  transitions.  The  answer  is  that  those  which  do  not  conduct, 
do  not,  at  the  critical  temperature,  become  displacive.  From  the  point  of 
view  of  that  most  basic  of  theories:  fluctuation-dissipation  theory,  a 
simple  answer  is  that  when  the  vibrational  wave  vector  becomes  larger  than 
the  correlation  length,  displacement  occurs.  It  is  as  if  the  vibrational 
cup  runneth  over  (to  become  a  pump  for  ezcitonic  coupling).  For  (CH)^  type 
of  polymers,  this  can  only  occur  when  the  energy  of  the  normal  mode  of  the 
vibrational  levels  within  the  electronic  potential  well  becomes  a  magnitude 
that  phonon-exci tonic  lateral  coupling  is  possible.  Parametric  excitation 
then  results,  which  is  describable  after  translation,  by  Eqs.  (19)  and  (24). 
The  vibrational  normal  mode  acts  as  a  pump  to  the  exciton  as  signal  and  the 
vibrational  normal  mode  of  the  neighboring  well  as  idler.  In  the  pendulum 
analogy,  the  vibrational  energy  levels  correspond  to  applied  torque,  in  the 
long  Josephson  junction,  to  the  applied  current.  The  excitonic  coupling 
corresponds  in  the  pendulum  analogy  to  the  torsion-bar  spring  constant,  and 
in  the  long  Josephson  junction  to  q/2jtL.  The  analysis  of  the  previous 
section  thus  extends  the  pendulum  analogy  to  a  consideration  of  the 
resultant  effect  of  vibrational  fluctuations  in  the  applied  torque. 

The  dynamics  of  fluctuations  in  a  critical  region  is  dominated  by  the 
role  played  by  exchange  of  energy  between  the  different  hydrodynamic  (long- 
wavelength)  motions  and  the  thermal  reservoir  (short-wavelength  motions). 
For  a  small  deviation  from  equilibrium,  the  rate  of  change  of  is 


proportional  to: 
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where  f  is  i  kinetic  coefficient.  As  before,  the  ordering  field  is  of  the 


F  «  /[j  a4>*  +  J  |  C(^4>)»  -  h4>]dnx,  (47) 

and  thus  the  equation  of  motion  (46)  has  the  form  with  c=l  and  a=r 

(Patashinskii  and  Pokrovskii,  1979;  Eilenberger,  1981): 

f  If  "  Vl<j)  +  T(f)  +  =  h  (48) 

Eq.  (48)  describes  the  situation  depicted  in  Fig.  9D  of  section  VI  for  which 

<J>  is  conserved,  E  is  conserved  and  F  =  -  A  V, 

o  o 

For  the  case  of  A  (4>  and  E  are  not  conserved  and  is  a  constant),  Eq. 


(46)  is: 


(t  $  —  V *4>  +  4g  $*)  +  h 

o  o 


For  the  case  of  B  (<j>  is  conserved,  E  is  not  conserved  and  r  =  -  A  V) 

o  o 


Eq.  (46)  is: 


&  +  $  t  .  o 

3 1  ‘  J 


where  the  flux  j  is  defined:  j  *  r  ^(f  -(6H/S4>)  and  f  (x,t)  is  f< 

o  ext  ext 

externally  applied. 

For  the  case  C  ($  is  not  conserved,  E  is  conserved  and  T  is  a 

0 

constant),  Eq.  (47)  is  (Patashinskii  and  Pokrovskii,  1979): 

H  “  To^*  +  2  (^)a  +  8o<i)4  +  +  2  °ol  e*-|dIlx' 

where  two  independent  fields  $  and  e  are  postulated  obeying: 

It  ♦«*.*>. 

It  ■-  4o  7'  S  +  T  <I-tK 


(52a) 


(52b) 


The  case  generalized  to  two  fields  <|>^  and  <j>2»  where  is  an  m-component 
field  and  $2  is  an  n-component  field  is: 

H  -  /ddx[  |  (t14>1*  +  T24>2»)  +  J  (V^)*  +  |  (V<t>2»)  + 

8,(^1*)*  +  go  (4>~* ) 2  +  g,- 


(53) 
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Such  a  system  has  two  phase  transitions  for  which  the  renormalization  group 
equations  in  four  space  are  (Patashinskii  and  Pokrovskii,  1979): 
if1  =  -4(m+8)gl*  -  ng12», 

if2  =  -4(n+8)g2»  -  mg12>,  (54) 

Jfl2  =-412((m+2)gl  +  (n+2)g2)  -8.^*. 

This  analysis  also  gives  a  relation  to  the  susceptibility  .  Returning  to 
Eqs.  (46)  and  (47),  we  obtain  the  form: 

$  =  -  T(a4>  +  b4> 3  -  CV  *4>  -  h)  .  (55) 

If  4>  is  close  to  the  equilibrium  value  <J>  ,  this  equation  can  be  linearized 
(<|>-<J>  +  64>) : 

Y  -  cV*6<t>  +  ^  =0.  (56) 

Changing  to  the  Fourier  components: 

d 

5<f.(q)  =  /6< Mx)e“1,,ld  x,  (57) 

Ve  see  that  their  relaxation  time  depends  on  q: 

t-Mq)  =  r/X(q).  (58) 

c 

and  the  susceptibility  for  a  given  q  is  determined  by  the  relation 

(Patashinskii  and  Pokrovskii,  1979): 

X~x( <l)  “  X“x(0)  +  cq*  •  (59) 

Could,  then,  any  arbitrary  polymer  system  exhibit  soliton  propagation? 

Aaly  if  it  has  the  requisite  hydrodynamic  vibrational  and  relaxational  modes 

which  occurs  when  the  wavenumber,  q,  and  frequencies,  a),  of  the  hydrodynamic 

vibrational  and  relaxational  motions  are: 

q  <<  1/r  ,  (60) 

c 

03  >>  1/t  , 

C 

where  r  is  the  static  correlation  length  and  t  is  the  relaxation  time, 
c  c 

If,  then,  we  divide  a  trans-polyacetyl ene  chain  into  cells  with  linear 

dimensions  dx  <<  r  ,  and  stipulate  that  in  each  of  these  cells  the  total 
c 


-35- 


NADC-83074-60 


order  parameter  varies  with  time,  the  characteristic  time  t(dx)  has  a  well 
defined  scaling  dimension  dt: 

t(dx)  ~  dx  dt,  (61a) 

and  the  inverse  co(dx)  =  t-x(dx)  is  the  characteristic  frequency  of  a 

fluctuation  of  size  dx  (i.e.  of  the  size  of  one  monomer  of  (CH)  ,  For  a 

x 

Q.  1/ *  i.e.,  not  the  hydrodynamic  but  the  f luctuat ional  mode,  the 

c 

dispersion  of  the  f luctuational  mode  is: 

-dt 

w  ~<L  (61b) 

Location  of  the  hydrodynamic  (both  vibrational  and  relaxational  modes)  and 


the  flue tuat ional  mode  is  shown  below  (after  Patashinskii  and  Pokrovskii, 


1979) : 


q 


i 

FLUCTUAT  IONAL 
MODE 


VIBRATIONAL' 
V ;  MODES/, 

</zzzs< 


0,0 

. 

THERMAL  CONDUCTION 

'(RELAXATIONAL  MODES) 


Fig .  15 

The  hydrodynamic  modes,  the  vibrational  and  relaxational,  are  associated 
with  the  transport  of  additive  conserved  quantities  such  as  spin  and 
electric  charge.  In  the  case  of  the  vibrational  mode,  there  is  almost  no 
dissipation,  so  case  D  above,  belongs  to  this  category  as  does  a  single  line 
of  trans-polyacetyl ene .  In  the  case  of  the  relaxational  mode  fluxes  are 
proportional  to  the  gradients  of  hydrodynamic  quantities,  and  cases  A,  B, 
and  C  above,  belong  to  this  category. 
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IX.  CHANGING  AND  INCOMPLETE  DEGENERACY 
a.  Introduction 

Trans-polyacetylene,  (CH)^,  is  electrically  conducting  and  the 
underlying  mechanism  mediating  this  conduction  is  thought  to  be  due  to  the 
presence  of  dynamic  force  relations  describable  by  the  mathematical  entity 
of  the  soliton  (Heeger  and  MacDiarmid,  1981).  There  are  a  number  of 
current-voltage  relations  reported  for  (CH)^  which  show  a  quadratic 
dependence  of  current  on  the  voltage  (Shirakawa  et  al.,  1978;  Park  et  al., 
1980;  Epstein  et  al.,  1980;  Roland  et  al.,  1980;  Gould  et  al.,  1981;  Hoses 
et  al.,  1981;  Chiang  and  Franklin,  1981;  Yamamoto  et  al.,  1980)  (cf.  Figure 
1).  With  the  dopant  concentration  y  >  10— a  the  EPR  signal  is  significantly 
reduced  from  that  of  the  undoped  (and  electrically  nonconducting)  trans- 
(CH)^;  and  for  y  >  2  x  10-1  the  signal  vanishes  (Ikehata  et  al.,  1980; 

Chien  et  al.,  1982),  indicating  that  the  dopant  binds  to  free  electrons 
providing  the  conditions  for  electron-free,  charged  excitation  (soliton) 
propagation.  This  doping  is,  however,  not  expected  to  be  uniform  throughout 
any  (CH)^  chain. 

There  have  been  a  considerable  number  of  theoretical  treatments  of 
soliton  propagation  in  (CH)^  (Su  et  al.,  1979,  1980;  Eivelson,  1981a, b).  A 
different  point  of  view  is  considered  here  permitted  by  the  NMR  observations 
of  Nechtschein  et  al  (1980)  on  changes  in  the  proton  relaxation  rate  upon 
doping  of  (CH)^  and  by  Chen  et  al.,  (1978)  that  doping  adds  to  the  double 
bond  and  at  the  highest  concentration  results  in  the  formation  of  a 
substituted  polyethylene.  These  observations  suggest  that  the  doping 
process  results  in  a  competition  between  the  dopant  and  the  hydrogen  (state 
of  resonance)  for  a  carbon  bond. 
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Fig.  16.  Current-voltage  characteristics  of  trans-(CH) 

at  various  temperatures.  From  (Yamomoto  et  al., 

1980) . 

As  electrical  conductivity  in  (CH)^  is  dependent  upon  this  doping, 
considered  nonuniform,  the  possible  reasons  are  examined  here  for  the 
quadratic  current-voltage  relations  in  trans-CCH)^  using  a  soliton  charge 
carrying  mechanism  based  on  the  enharmonic  oscillator  model  (Chen  et  al., 
1978). 

b.  Anharmonic  Oscillator  Model 

The  anharmonic  oscillator  model  is  based  on  the  equation  of  motion 
for  a  monomer  of  the  general  form  (Bruce,  1978;  Bruce  and  Cowley,  1980): 

i  -  x*  +  gX  +  C  «  0,  (62a) 

where  x  is  the  B  position  of  the  carbon  bond  and  <t>  is  a  configurational 
potential  (Figure  2),  and  for  coupling  to  2d  immediate  neighbors: 

”  I  a|  */4B  is  the  depth  of  the  local  potential,  (62b) 

c|  A|  /Bd  is  the  representative  bond  energy. 


(62c) 
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defining  a  constant  h®E^/Eg.  (62d) 

For  h>>l  only  localized  transitions  occur,  but  for  h<<l  a  displacive 
transition  (Brnce  and  Cowley,  1980;  Murata,  1975;  Anbrey,  1975,  1976; 


X 

Fig.  17 

McMillan,  1976;  Hohenberg  and  Halperin,  1977;  Halperin,  Hohenberg  and  Ma, 
1974;  Halperin,  Hohenberg  and  Ma,  1976)  results  in  an  incommensurate  phase 
domain  wall  diffusion  (soliton  propagation). 

To  describe  a  propagating  transition  in  a  chain  of  monomers  of  the 
type  depicted  in  Figure  17,  the  transition  becomes  displacive  (Schneider  and 
Stoll,  1973,  1975,  1976,  1978)  (h<<l),  resulting  in  domain  wall  diffusion. 
The  potential  energy  of  the  Hamiltonian  description  of  the  process  is 
(Bruce,  1978;  Bruce  and  Cowley,  1980): 

V-  §  ***<u(£)  -  uU*)»)  +£$(£)),  (62e) 

~  X.  £ 

Taking  into  consideration  that  the  doping  process  results  in  the 
previously  mentioned  competition  with  the  hydrogen  (state  of  resonance)  for 
a  carbon  bond  (Chen  et  al.,  1978),  then  ideally  the  doping  results  in  two 
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types  of  configuration,  A  and  B  (Figure  18).  In  the  Figure  18A  doping 
condition  the  bond  transfer  from  left 
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to  right  (effecting  a  transfer  from  one  isomerized  state  to  the  other)  is 
facilitated,  as  the  dopant  provides  a  stronger  electron  acceptance  in  the 
isomerized  state  to  be  adopted.  In  the  Figure  18B  doping  condition,  on  the 
other  hand,  the  bond  transfer  from  left  to  right  is  less  likely,  as  the 
dopant  provides  a  stronger  electron  acceptance  in  the  isomerized  state 
already  existing.  The  actual  result  of  doping  is  more  probably  the  nonideal 
forms  shown  in  Figure  19,  where  the  angle  X  is  a  pictorial  representation 
for  the  degree  of  resonance  interaction  of  bond  1  with  bond  2  and  the  depth 
of  the  configurational  double  potential  well  (Figure  20). 
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It  is  a  consequence  of  the  enharmonic  oscillator  model  that  any  of 
the  joined  monomer  segments  shown  in  Figure  19  may  be  represented  as  a 
trajectory  over  the  a,b  surface  representation  of  the  monomer  element 
dynamics  described  by  Eq.  (62a).  In  Figure  21  is  shown  an  arbitrary 
semicircular  trajectory;  and  Figure  22  depicts  the  possible  nonuniform 
doping  of  (CH)^  and  monomer  configurational  potential  wells  for  a  complete 
circular  trajectory.  The  particular  trajectory  is,  of  course,  arbitrary, 
but  for  purposes  of  exposition  the  circular  trajectory  will  be  used  here. 

In  the  case  of  trajectory  l-*-8,  at  each  of  the  8  monomer  sites  and  for  an 
increasing  applied  current  to  the  (CH)^  chain,  the  measured  voltage,  if  it 
could  be  measured,  would  show  an  initial  increase  to  V  and  then  greatly 
diminishing  increases  after  that  value.  Conversely,  for  an  increasing 
applied  voltage  to  the  (CH)^  chain,  the  measured  current,  if  it  could  be 
measured,  would  show  an  initial  slow  increase  until  a  AI  value  is  reached 
and  then  increase  greatly. 

U. 


H  H  H  D  H  H 

®  ®  <D  ®  ® 


Fig.  21 
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With  constant  current  applied  and  the  resistance  and  voltage 
measured  at  each  monomer  location,  if  snch  could  be  measured,  the  following 
relation  wonld  then  be  obtained  (Figure  23): 


CONSTANT  APPLIED  CURRENT 


Fig.  23.  Logarithmic  representation  of  AV,  1/AR  and  constant  applied  AI 
relations  of  the  trajectory  of  Figure  22.  The  numbers  refer  to 
the  locations  in  Figures  21  and  22. 

For  the  particular  nonuniformly  doped  chain  considered,  the  monomer  current 
voltage  relations  at  each  monomer  site  would  be  (Figure  24): 


2  3  1  4  8  75  6 


Fig.  24.  Current  singularities  in  the  nonuniformly  doped  (CH)^  chain  for 
the  locations  in  Figures  21  and  22. 
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There  is  little  hope  of  measuring  such  monomer  current-voltage 
relations  presently.  However,  the  task  to  relate  the  sol iton-mediated 
conductivity  at  each  monomer  site  to  the  conductivity  which  may  be  measured 
for  a  collection  of  such  chains  may  be  approached  analytically.  For 
example,  ( CH)^  chains  consisting  of  monomers  with  less  resistance,  AR, 
should  exhibit  more  solitons,  at  equal  applied  current,  than  chains  with 
more  such  resistance. 

It  is  demonstrated  in  the  next  section  that  relating  the  current- 
voltage  quadratic  relation  problem  to  the  d.c.  current  singularities  of  the 
long  Josephson  junction  is  straightforward,  resulting  in  an  explanation  for 
the  current-voltage  relations  of  a  collection  of  nonuniformly  doped  (CH)^ 
chains,  as  both  phenomena  are  based  on  the  same  soliton  dynamics, 
c.  Long  Josephson  Junction  Analogy 

The  problem  of  predicting  the  average  voltage  for  a  given  applied 
current  in  a  nonuniformly  doped  (CH)^  chain  is  directly  related  to  the 
explanation  of  the  d.c.  current  singularities  of  the  long  Josephson  junction 
(Chen  et  al.,  1971;  Parmentier,  1978;  Parmentier  and  Costabile,  1978; 
Costabile  et  al.,  1978;  Costabile  et  al.,  1980;  Costabile  and  Parmentier, 
1980)  (Figure  25).  The  fluxon  on  the  long  Josephson  junction  is 
mathematically  equivalent  to  a  soliton  and  shares  with  the  (CH)^  soliton  a 
common  analytical  mechanics. 
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Fig.  25.  Current-voltage  characteristic  of  a  very  long  Josephson  junction 
showing  zero— field  current  steps  in  the  presence  of  an  external 
magnetic  field  applied  transversely  across  the  junction.  From 
Costabile  et  al.,  1980. 

The  sine— Gordon  equation,  without  dissipative  effects  for  which 


the  soliton  is  a  solution,  is: 


3»<b  _  a»d) 


sin<j>. 


(62e) 


3x*  at* 

where  $  is  the  phase  difference  between  macroscopic  quantum  wave  functions 
that  characterize  two  superconductors  and  is  directly  analogous  to  the  4>  of 
the  enharmonic  oscillator  model.  The  soliton  (fluxon)  solution  to  Equ. 


(62e)  is  (Lamb,  1971): 

<p  ■  4  tan-*(f(x)g(t)) , 


(62f) 
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and  f  and  g  are  Jacobian  elliptic  functions  defined  (Chen  et  al.,  1971; 

Parmentier,  1978;  Parmentier  and  Costabile,  1978;  Costabile  et  al.,  1978): 
df  * 


dz 
da  1 
dt 


=  Af4  +  (1  +  B) f *  -  C,  (62g)  (i) 

=  Cg4  +  Bg*  -  A,  (62g)  (ii) 

where  A,  B  and  C  are  arbitrary  constants.  Taking  a  =  (1+B)/A,  b  =  C/A  for 
Equ.  ( i)  and  a  =  B/C,  b  =  A/C  for  (ii),  and  differentiating,  the  form  of 
Equ. (62a)  is  obtained. 

For  the  case  of  most  interest,  that  with  dissipative  effects,  we  have: 

(62h) 

where  a  and  £  are  constants.  This  form  can  be  solved  by  assuming  an  z  and  t 
dependence  of  <p,  e.g.: 

<p  =  <f>  (z  -  ut)  =  <£(£),  (62i) 

where  a  is  a  constant  propagation  velocity. 

With  o:  and  £  related  to  an  external  energy  source,  Equ.  (62g)  is: 

(62  j) 


dT*  "  ~  a  11  +  cos^  ft  =  sin*' 


124,  _  124  -  r  (1  +  cos<|>)  ||^  =  sin  -  y, 


dx* 

where  F  =  g  $  /2nC, 
o  o 

£  = 

1  o 


3t* 


Y 


vv 


a  -  g  v  ($  /2nl  C)  V* 
o  o  o 

and  I_  is  a  uniformly  distributed  current  bias  for  the  long  Josepbson 

D 

junction  and  a  current  bias  for  the  (CH)^  chain. 

With  this  external  source  dependence,  Equ.  (62j)  has  the  desired 
quadratic  current-voltage  characteristic  with  normalized  voltage 
(Parmentier,  1978): 

V  IV 

(62k) 


<*?>*/■  da  1 1  c0>;  H. 


where  is  an  arbitrary  constant; 


(x  -  ut) ; 


and  dn  [;]  represents  a  Jacobian  elliptic  function  (Lamb,  1971)  of  modulus 
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k,  0<k<l; 

k=  <^77-  >v»;  (621) 

Y  -  2^ u*  [ ( 1— u* ) 1  +  4r*u*]V*.  (62m) 

o 

Because 

V  (  =  2nK(k),  (62n) 

k  21  u* 

where  i  is  the  length  of  the  long  Josephson  junction  (or  (CH)^  chain; 

K(k)  is  the  complete  elliptic  integral  of  the  first  kind;  and  n  is  the 
number  of  solitons  (fluzons  contained  on  it),  the  average  value  of  voltage, 
v,  is: 


<v>  =  nu/£. 


(62o) 


As  the  number  of  solitons,  n,  will  vary  from  monomer  to  monomer  at  constant 
applied  current,  depending  on  the  doping  conditions,  the  average  <n>  is 
substituted  for  n  and: 

<v>  *  <n>u/£.  (62p) 

The  current-voltage  relations  calculated  with  Equ.  (62p)  are  shown  in 
Figure  26.  The  similarity  between  the  relations  shown  in  Figure  26,  the 


Fig.  26.  Current-voltage  characteristics  for  n  fluzon  oscillations.  From 
Parmentier,  1978.  Fluzon  oscillations  move  to  one  end  of  a  long 
Josephson  junction  are  reflected  as  an  antifluzon  and  move  to  the 
other  end,  etc.  (Fulton  and  Dynes,  1973).  Such  oscillations  can 
adsorb  power  from  an  ezternal  source  of  bias  current.  Costabile 
et  al.  (1978)  have  also  made  similar  calculations  using  an 
approzimat ion  scheme  rather  than  an  ezact  solution  with  periodic 
boundary  conditions. 
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putative  relations  of  Figure  24  and  the  experimentally  obtained  relations  of 
Figure  16,  indicate  that  a  possible  explanation  of  the  nmlinearity  of 
Figure  16  is  obtained. 

X.  THE  IMPORTANCE  OF  SYMMETRY  BREAKING 

In  previous  sections,  the  precise  underlying  dynamics  of  each  link  in 
the  polyacetylene  chain  was  addressed.  The  behavior  of  such  a  one¬ 
dimensional  conductor  in  T-junction  form  was  described  (section  I-III)  and 
the  effect  of  uneven  doping  and  the  application  of  Ohm's  and  Eirchoff's  laws 
to  these  systems  predicted  (sections  TV-IX) . 

In  this  section,  two  facets  of  the  one-dimensional  organic  conductor 
problem  are  addressed.  The  first  is  a  description  of  how  the  conditions  for 
the  underlying  dynamics  of  soliton  propagation  may  arise,  and  we  show  that 
these  conditions  are  caused  by  broken  symmetry,  placing  the  problem  squarely 
in  the  frame  of  modern  critical  phenomena  theory. 

The  description  of  this  first  facet  leads  to  a  description  of  the 
second,  which  is  the  relation  of  our  present  description  of  every  link  in 
the  polyacetylene  chain  to  the  more  manageable  Hubbard  model  (Hubbard,  1963, 
1964a, b)  and  its  developed  form  (Economou,  White  and  DeMarco,  1978:  White 
and  Economou,  1978:  DeMarco,  Economou  and  White,  1978).  In  achieving  a 
description  of  the  second  facet  of  the  underlying  dynamics,  the  problem  is 
placed  in  the  frame  of  current  correlation  techniques  for  describing 
conduction  and  magnetism  based  on  the  double-time  Green  function  method  of 
Zubarov  (1960)  (section  XI  below).  In  effect,  this  reduces  the  nonlinear 
many-body  problem  addressed  in  sections  I-IX  to  a  one-body  one. 

In  reducing  the  many-body  problem  to  one-body  (Hubbard)  form,  we  see 
advantages  and  disadvantages.  One  advantage  is  the  obtainment  of  a  more 
manageable  Hamiltonian  and  in  relating  this  development  to  mainstream 
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theories  of  conduction.  A  disadvantage,  however,  is  that  the  underlying 
dynamics  is  introduced  somewhat  artifically  to  the  essentially  static  one- 
body  description  with  a  resulting  decrease  in  physical  insight.  Thus,  we 
recommend  that  only  the  sections  I-IX  and  present,  taken  together,  achieve 
an  in  depth  description  of  the  underlying  dynamics. 

The  following  is  a  general  description  of  symmetry  breaking  at  both  the 
micro  and  macrophysical  levels  beginning  with  a  description  of  polymer 
monomer  described  by  the  Schrodinger-like  equation  (20)  of  section  TV: 

y  -  f*  +  af  +  b  =  0,  (62b) 

o  o 

where  f,  may  be  taken  as  the  midfrequency  of  a  signal  equivalent  to  electron 
momentum  in  the  analysis  of  states.  The  procedure  is  to  demonstrate  that  a 
system  described  by  an  equation  similar  to  that  of  Eq.  (62)  not  only 
exhibits  two  potential  minima  (isomerization)  when  the  variable,  a,  is 
negative  (which  is  known),  but  that,  more  importantly,  a  negative  variable, 
a,  results  in  symmetry  breaking. 

In  order  to  achieve  this  aim,  some  conclusions  are  first  summarized 
before  (Barrett,  1977,  1978,  1979a, b. ,1980)  pertaining  to  any  system 
described  by  an  equation  similar  to  Eq.  (62). 

a.  Although  never  stated  as  a  logical  consequence  of  any  deeper 
theoretical  condition,  the  fundamental  conditions  for  the  valid  application 
of  quantum  mechanics  are  the  commutation  relations  (Dirac,  1967:  Heitler, 
1954:  Schiff,  1973:  Louisell,  1973): 
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where  the  operators  q  and  p  are  canonical  coordinates  and  momenta  and  r  and 
s  refer  to  degrees  of  freedom,  r  =  s.  With  one  degree  of  freedom,  the 
commutation  relation  is: 

qp  -  pq  =  iA.  (64) 

b.  The  supposition  of  this  commutation  relation  enables,  using  the 
Schwartz  inequality,  the  following  uncertainty  condition  to  be  described 
involving  the  numbers,  Ap  and  Aq,  defined  as  standard  deviations  about  p 
and  q: 

Ap  Aq  =  1/2A.  (65) 

Because,  just  as  readily,  the  following  comutation  relation  might  be 
supposed  involving  the  operators,  Ap  and  Aq: 

AqAp  “  ApAq  =  ik,  (66) 

then  we  arrive  by  similar  reasoning  at  the  uncertainty  relation  for  the 
numbers,  p  and  q  (Barrett,  1977): 

p-q  =  Ulk,  (67) 

thereby  defining  the  quantum  condition  as  a  four-parameter  condition  in 
which  the  concept  of  number  or  operator  is  considered  relative. 

c.  Because  of  logic  used  to  solve  the  Schrodinger  equation  and  to 
define  the  energy  eigenfunctions  of  states  of  the  harmonic  oscillator  could 
as  easily,  and  in  the  long  run  more  appropriately,  define  signals  produced 
by  a  harmonic  oscillator  (i.e.,  the  uncertainty  relations  might  be  derived 
as  a  condition  for  the  solution  of  a  harmonic  oscillator  -  this  is  admitted 
as  possible  in  some  texts,  e.g.,  (Schiff,  1973,  p.  73),  we  have  then  the 
possibility  of  more  general  relations. 

Because:  (i)  In  the  root-mean-square  deviation  from  the  mean,  or 
expectation  value,  (Ap)1  =  <(p-<p>)2>  =  <p2>  -  <p>2, 

(Aq) 2  =  <  ( q- < q> )  * >  =  <q2>  -  <q>2. 
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and  (ii)  as  <p>  *  <q>  *  0  for  any  harmonic  oscillator  wave 
function,  then  (i)  permits  ns  to  write: 

(Ap)*  =  <p>*, 

( Aq) *  -  <  q> 1 . 

(iii)  If  flnctnations  are  now  introduced,  not  to  p  and  q  in  the 
sense  of  defining  Ap  and  Aq,  bnt  to  Ap  and  Aq  in  the  sense  of  defining  p  and 
q,  then,  in  analogy  to  (i): 

(p) *  -  <Ap  -  <Ap> ) * >  *  < Ap* >  -  <Ap>», 

(q)  *  =  <Aq  -  <Aq> ) * >  *  <Aq*>  -  <Aq>», 

i.e.,  if  p  and  q  are  thought  of  as  types  of  carrier  signals  and  Ap  and  Aq  as 
types  of  modulating  envelopes  of  those  signals,  then  (i)  above  describes 
fluctuations  in  the  carrier  signals  and  here  in  (iii),  fluctuations  in  the 
modulating  envelopes  are  described. 

(iv)  Thus,  in  analogy  to  (ii):  <Ap>  =  <Aq>  =  0  for  any  modulated 
harmonic  oscillator,  and  (iii)  permits  us  to  write: 

(p) »  *  <Ap>*, 

(q) »  *  <  Aq>  * . 

Ve  then  have  the  two  uncertainty  products  for  the  modulated  harmonic 
oscillator: 

Ap.Aq  =  l/2(2n  +  1)K:  p*q  =  l/2(2n  +  1)H,  n  =  0,  1,  2,  ...  (68) 

The  physical  significance  of  this  viewpoint  is  that  an  uncertainty 
relation  may  arise  either  from  inherent  fluctuations  in  the  oscillator  or  in 
fluctuations  in  the  external  modulation  of  the  oscillator. 

d.  With  an  emphasis  on  the  signal-like  character  of  the  uncertainty 
products,  which  in  acoustics  can  be  used  as  signals,  it  is  easier  to  treat 
fQ,  or  signal  midfrequency,  instead  of  p,  Af,  or  signal  bandwidth,  instead 
of  Ap,  t^,  or  signal  midperiod,  instead  of  q,  and  At,  or  signal  duration, 
instead  of  Aq. 
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This  permits  the  following  relation  for  qnantal  signals  in  a  way 
identical  to  the  uncertainty  product  derivation  of  c: 


f  •  t  =  l/2c(2n  +  1),  n  =  0,  1,  2,  ....  (69a) 

o  o 

Af*  At  =  l/2c(2n  +  1),  n  =  0,  1,  2,  ....  (69b) 

Af  =  f  /f,  At  =  4t  ,  for  n  =  0,  (69c) 

o  o 

Af*At/f  *t  =1,  (69d) 

o  o 


i.e.,  a  constant  magnetic  field  influence  is  assumed  (constant  modulation), 
n  is  the  angular  momentum  quantum  number,  c  =  Ji  for  wave  optics  or  reception 
by  a  radially  circulating  receiver,  and  c  =  1  for  wave  acoustics  o  a 
linearly  acting  receiver. 

e.  The  wave  packet  or  qnantal  signal  for  light  and  sound  with  n  -  0 


is: 

s(t)  =  exp[-(l/c)(t  -  t  )2]exp[i2nf<i>t]  in  the  time  domain,  (70a) 

o 

S(f)  =  exp[-c(n)a(f  -  fQ)a]exp[-i2ntof]  in  the  frequency  domain.  (7v.b) 
In  general,  the  signals  are: 

s(t)  =  D  (t)cos (180t/%(2n  +  1))°,  n  =  0,  1,  2,  ...»  (71) 

n 

where  the  D  (t)  are  modified  Weber-Hermite  polynominals : 
n 

D  (t)  «  exp[(-l/2c)t»]H  (t),  n  -  0,  1,  2,  ...  (72) 

n  n 

and  the  HQ(t)  are  Hermite  polynomials: 

H  (t)  =  (-l)nexp[t*]3nexp[-t»]/3tn,  n  -  0,  1,  2,  ...  (73) 

n 


The  equations  expressed  in  Eqs.  (69)  are  the  solutions  to  the  wave 
amplitude  equation  for  the  harmonic  oscillator  system  (i.e.,  to 
Schrodinger ' s  equation).  A  formally  similar  equation  to  the  Schrodinger 
equation  relates  all  four,  but  only  those  four,  signal  parameters  has  been 
shown  to  be  (Barrett,  1977): 

f  3  -  (-4Af»At/t  )f  +  (-16Afa/t  )  =  0,  (74) 

O  O  O  0 
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This,  in  turn,  is  example  of  the  well-studied  equation  for  a  harmonic 
oscillator: 

$  =  f  3  +  af  +  b  =  0,  (75) 

o  o 

or  the  differential  of: 

<p  =  1/ 4f  4  +  l/2af  *  +  bf,  (76) 

o  o 

where  <J>  is  a  potential  function. 

f.  In  Figs.  28  and  29  curves  of  <{>  vs.  f  are  plotted,  demonstrating 
that,  within  a  certain  range  of  a  =  4Af*At/tQ  and  b  =  -16Afa/tQ  values,  two 
potential  minima  are  available,  and  the  quantum  effect  is  inseparable  from 
the  existence  of  an  area  of  instability  defined  by  the  presence  of  more  than 
one  potential  minimum. 
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Here,  the  intention  is  to  demonstrate  that  these  two  potential 
minima  define  the  location  of  tvo  fired  points:  the  Gaussian  and  a 
nontrivial  fixed  point,  where  a  fixed  point  is  defined  in  the  terms  of 
renormalization  group  theory  (Pfeuty  and  Toulouse,  1977:  Wilson,  1971)  as  a 
point  in  parameter  space  vhich  by  itself  constitutes  its  own  entire 
trajectory.  To  such  a  point  there  corresponds  a  state  of  the  system 
invariant  under  the  operations  of  the  renormalization  group  and  which  has  a 
correlation  length  either  infinite  or  zero  (Pfeuty  and  Toulouse,  p.  13). 

( i) .  Argument 

In  the  vicinity  of  a  transition  point  between  a  stable  region 
(defined  by  one  fixed  point,  the  Gaussian)  and  an  unstable  region  (defined, 
as  above,  by  two  fixed  points,  one  of  which  is  the  Gaussian),  there  are 
fluctuations  which  govern  the  behavior  of  the  system.  Close  to  the  critical 
point  quantities  may  be  observed  which  obey  power  laws  with  exponents  that 
are  not  integers.  For  example,  just  below  Tc,  the  critical  temperature,  the 

O 

order  parameter  (f  )  begins  by  falling  proportionately  to  (T  -  T)  ,  while 
o  c 

just  above  T  the  specific  heat,  for  instance,  is  proportional  to  (T  -  T  )  °, 
c  c 

-V 

the  correlation  length  to  (T  -  T  )  ,  and  the  susceptibility,  in  the  case  of  a 

c 

magnetic  system,  to  (T  -  T  )~ ^.  The  exponents  a,  P>  v,  y,  and  others  are 

c 

defined  in  the  vicinity  of  a  given  critical  point,  and  are  called  critical 
exponents.  These  critical  exponents  obey  some  simple  relations  called  scaling 
laws,  which  have  a  degree  of  universality.  In  the  case  of  a  magnetic  system, 
the  scaling  laws  are  satisfied  automatically  if  the  singular  part  of  the  free 
energy  G(T  -  Tc<  H)  and  the  correlation  function  T(T  -  Tc»  R  )  are  homogeneous 
functions  of  their  arguments. 

In  the  case  of  a  magnetic  system,  for  the  equation  of  state 
there  exists  a  law  of  corresponding  states  interrelating  the  order  parameter. 
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M,  it*  conjugate  field,  H,  and  the  temperature,  (T  -  T  =  r)  by  an  equation  of 

c 

the  form  (Stanley,  1980): 

H/M3  -  h(T/M1/P)  (77) 

In  the  case  of  the  quantal  system  defined  by  Equs.  (69),  this 
is: 

t  /f  5  =  lu  At/ f  1/P) .  (78) 

o  o  o 

With  8  =  1  and  substituting  c  =  1  for  h,  the  relations  expressed  by  Equs.  (69) 

and  (78)  give  a  value  for  0  of  log(f  )/log(4f  ),  e.g.  0.6  for  f  =  8.  In  the 

0  0  o 

case  of  the  relations  expressed  by  Equs.  (1),  0  =  log(p)/ log(4hp) . 

In  terms  of  a  modulated  harmonic  oscillator,  Equ.  (78) 

describes  the  relation  of  the  midperiod  of  oscillation  (fQ),  the  midperiod  of 

oscillation  (tQ),  and  the  temporal  window  of  modulation  (  t),  for  the  case  of 

the  minimum  quantum  (n  =  0,  Equ.  (69)  at  different  0s  (*  log(f  )/log(4f  ))  or 

o  o 

different  fQs.  A  similar  description  holds  for  p  and  Ap  for  the  case  of  the 
minimum  quantum  (n  =  0,  Equ.  (68))  and  different  p's. 

From  Widom's  scaling  law  (Pfenty  and  Toulouse,  1977,  p.  22): 

Y  **  0(6  -  1),  (79) 

we  have,  for  the  present  system,  y  ■  0.  From  Rushbrooke's  scaling  law: 

a  +  20  +  y  -  2,  (80) 

we  obtain  a  =  0.8  for  fQ  =  8,  or  u  =  2  -  21og(fQ)/ log(4fQ)  generally,  and  from 
Fisher's  scaling  law: 

Y  -  (2  -  n)v,  (81) 

either  v  =  0  or  n  ■  2.  As  Josephson's  scaling  law: 


vd  =2-o 


(82) 
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must  be  satisfied,  n  =  2.  With  d  =  4,  which  corresponds  to  a  -  0  (of  Equ. 

(62):  f  1  +  af  +  b  =  0),  Equ.  (82)  gives  v  =  0.3  for  f  =  8,  or  v= 
o  o  o 

l/21og(fo)/log(4fo)  generally.  Within  the  area  of  instability,  then,  a  <  0, 

<  4  and  v)  l/21og(fo/ log(4fQ) ,  for  the  qnantal  system  considered  with  the 
correlation  length,  £,  defined: 

<  -  imt;  <83> 

Josephson's  law,  vd  =  2  -  a,  links  the  space  dimensionality, 
d,  the  exponent,  a,  referring  to  a  thermodynamic  potential  and  the  exponent  v 
referring  to  a  correlation  function.  As  d  decreases  through  the  value  d  =  4, 
one  moves  from  a  regime  (d  >  4,  a  >  0),  where  Josephson's  law  fails,  to 
another  regime  (d  <  4,  a  <  0),  where  it  holds.  The  mathematical  basis  for 
this  lies  in  the  simple  identity,  -a  =  4  -  d  (where  a  is  defined  in  Equ.  (62) 
or,  as  is  generally  given,  in  the  expansion  parameter  e  =  4-d  (Hohenberg  and 
Halperin,  1977),  i.e.,  e  =  -a. 

The  relation  of  the  space  dimensionality  d  to  the  parameter  a 
of  Equ.  (74)  is  seen  clearly  in  Fig.  30  and  in  recent  developments  in 
renormalization  group  methods  (as  reviewed  in  (Hohenberg  and  Halperin,  1977). 
These  methods  include  1)  reduction  of  the  number  of  degrees  of  freedom,  2) 
restoration  of  the  spatial  density  of  degrees  of  freedom  to  its  initial  value 
and  3)  renormalization  of  spin  magnitude,  if  appropriate  to  do  so.  Referring 
to  Fig.  31,  a  comparison  of  two  systems,  1  and  2,  can  be  obtained  if  a,  is 
renormalized  to  equal  a.  (for  the  systems  f  +  a.f  +  b,,  i  *  1,  2).  The  two 
systems  would  then  be  located  at  the  same  position  along  the  a  axis  and 
comparisons  of  the  two  systems  could  be  made  along  the  fQ  and  b  axes.  This 
is,  then,  the  essence  of  the  renormalization  approach. 

If  elementary  oarticle  behavior  is  considered,  then  the  particle 
only  exists  in  the  region  d  <  4,  a  <  0,  due  to  the  energy  minimization 
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requirement  for  the  solution  of  a  harmonic  oscillator  given  in  X(c).  To 
ask,  therefore,  whether  an  elementary  particle  could  exist  in  the  classical 
domain  (d  >  4,  a  >  0)  is  to  ask  whether  the  energy  minimization  requirements 
could  be  disregarded. 

A  formal  analogy  can  be  developed  for  another  example  using 

the  solution  to  the  percolation  problem  by  Fortuin  and  Kasteleyn  (1972).  This 

application  illustrates  the  generality  of  the  mathematical  system  generated  by 

broken  symmetry,  i.e.  the  dynamic  underlying  both  the  quantum  harmonic 

oscillator  and  the  percolation  problem  is  formally  similar.  The  percolation 

problem  considers  an  infinite  lattice  whose  sites  are,  at  random,  either 

allowed  for  forbidden,  with  a  probability,  p,  that  any  given  site  is  allowed. 

A  critical  value,  p  ,  is  defined  so  that  when  p  2  P  *  there  occurs  an 

c  c 

infinitely  large  cluster  of  allowed  sites  and  percolation  takes  place,  i.e., 

one  can  cross  the  lattice  by  going  successively  from  one  allowed  site  to  a 

neighboring  allowed  site.  In  the  present  instance,  and  also  in  the  case  of 

ferromagnetism,  we  let  (p  -  p.  )  =  (t  -  t  )  »  At.  Other  formal  mathematical 

»  o 

correspondences  are: 

a.  The  probability,  N,  that  an  allowed  site  belongs  to  an 

infinite  cluster  is  equivalent  to  the  frequency,  f  ,  in  the  present  case,  and 

o 

to  spontaneous  magnetization  in  the  ferromagnetic  case, 

b.  The  average  size  of  finite  clusters  is  equivalent  to  the 
susceptibility  in  the  ferromagnetic  case  and  will  be  shown  to  be  unity  in  the 
present  case,  and 

c.  The  average  number,  G,  of  clusters  is  equivalent  to  the 

free  energy  in  the  ferromagnetic  case  and  will  be  shown  to  have  a  value  in  the 


present  case  when  n  ^  0. 
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For  the  general  case  expressed  by  Eqs.  (7)  and  (9),  the 

average  number  of  signals  of  size  n  is: 

G(t,n)  =  ,  (84) 

and  the  average  number  of  signals  having  any  size  is: 

G(t)  -  Eg( t,n) .  (85) 

n 

Then,  the  average  size  of  finite  signals  (susceptibility)  is: 

X  (t)  =  En*G(t,n) ,  (86) 

n 

or: 

x(t)  ~  (Jt)Y,  (87) 

whence  ve  see  that  for  G  and  X  to  be  other  than  unity,  y  =  0  or  5  ^  0  or  n  ^ 

0. 

The  percolation  problem  solution  also  gives: 

f  -  AtP,  (88) 

o 

which  is  satisfied  by  Equ.  (69)  for  8  =  log(fo)/ log(2/fQ) . 

We  turn  now  to  the  role  of  broken  symmetry  in  the  transition 
from  classical  to  quantum  mechanics.  The  Ginzburg  criterion  describes  the 
existence  of  a  characteristic  dimensionality  below  which  fluctuations  are  so 
large  that,  in  the  vicinity  of  the  critical  point,  the  self-consistency  of  the 
classical  theory  of  mechanics  is  destroyed,  resulting  in  a  broken  symmetry  for 
certain  rotations  and  the  necessity  of  a  quantum  mechanical  formulation.  This 
amounts  to  defining  B,  the  coefficient  of  the  quartic  term  in  the  usual 
expression  for  free  energy  (F  =  A*  Ma  +  B*  M4  +  k(VM)*),  as: 

B  ~  (At)U-d)/2.  (89) 

On  the  other  hand,  we  have  previously  defined  the  area  of 
instability  (in  which  two  potential  minima  exist,  not  one  minimum)  as 
extending  from  b  =  +1.3  (cf.  Fig.  28),  where  b  is  defined  in  Equs.  (75)  and 

3 

(76)  (f  +  af  +  b  ■  0) .  In  the  normalized  units  of  Fig.  28,  At  is  0.147, 
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and  for  d  =  4,  5,  6,  7,  8,  . ...  we  have  B  =  1,  3,  7,  18.5,  46.3,  ... 

To  describe  the  Hamiltonian  of  a  system  capable  of  broken 


symmetry,  we  let  At  ~  iq  (more  generally,  rQ  ~  At((l  +  At  +  ...)  + 

y  / / ) 

(At)  /v  7 (1  +  At  +  ...)))  and  B  ~  uq.  The  system  Hamiltonian  is  then: 


H  =  |t  +1/2  / 


(r  -  k»)  f(k)  *  + 


k  <  1/2 


»»«  /  f(k1).f(ka)*f(k,)S(k1  +  k»  +  k,  +  k4 ) . 

ki»k».k,,k4  <  1/2 

where  a  is  a  unit  of  length  for  the  system,  or  lattice  spacing.  Any 
dilatation  of  this  on  it  of  length  requires  renormalization  of  the  Hamiltonian, 
and  the  parameters  r^  and  uq  renormalize  (Wilson,  1971)  subject  to  the 
differential  equations: 


— r-  *  2r  +  buo(l  -  r  ), 
dX/  o  o 


(91a) 


dT  =  8X\>  "  CV 


(91b) 


where  e  ■  4  -  d,  b  =»  16(n  +  2),  and  c  =  16 (n  +  8),  and  s  =  e  ,  for  the 

dilatation  a  •+  sa.  Equs .  (91),  set  to  zero,  yield  two  solutions.  The  first, 

(r  *  =  u  *  *  0),  corresponds  to  the  Gaussian  fixed  point,  and  the  second, 
o  o 

(r*=-!r  e:u  =~),  corresponds  to  the  nontrivial  fixed  point.  In  the 
o  2c  o  c 

neighborhood  of  each  fixed  point  r  =  r  *  +  8r  and  u  =  u  *  +  6u  permitting 

oooooo 

the  definition  of  a  matrix  A  for  the  Gaussian  fixed  point: 

8 

2,  b 
A  -  (  ) 

g  0,  e 


with  eigenvalues  y^  =  2  and  =  s ,  and  a  matrix  A^^  for  the  nontrivial  fixed 


points : 


b  b 

(2  -  c£  ),  b(l  +  2cE  ) 
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with  eigenvalues  y..  =  2  -  —  and  y,  =  -a .  For  the  simplest  system  des- 
1  c  z 

cribed  by  Equ.  (69),  n  =  o  and  y^  =  2  -  e/4  and  y^  =  -e,  in  the  case  of  the 
nontrivial  fixed  point. 

This  development  permits  an  important  analogy  between  wave 
mechanics  and  lattice  behavior  i.e.,  the  underlying  dynamics  of  these  two 
distinct  physical  phenomena  have  a  similar  mathematical  form.  The  quantum 
system  we  are  considering  is  based  on  an  expansion  in  terms  of  n,  the 
angular  momentum  quantum  number  (cf.  Equs.  (69),  (70),  (71),  and  (72)). 
Another  way  of  considering  this  description  is  to  identify  n  with  a  lattice 
symmetry  index  or  degree  of  isotropy  (Fisher,  1974).  In  this  alternative 
description,  we  commence  with  a  lattice  with  sites  x  populated  with  spins, 

S^,  and  we  suppose  the  spin  vector  has  n  components,  i.e.,  =  (S^^),  p  = 

0,  1,  2,  ...  n,  which  enter  equally  into  interactions.  The  case  n  =  0 
describes  the  statistics  of  a  self-avoiding  walk  or  polymer  chain  in  solution 
(de  Gennes,  1972),  n  =  1  described  the  uniaxial  or  Ising  model,  n  =  2,  the 
XT  or  planar  model,  and  n  =  3,  the  Heisenberg  model.  The  general  model  is 
the  eight-vertex  model  of  Baxter  (Baxter,  1971).  In  this  case,  there  are 
eight  possible  different  configurations  of  arrows  (spins)  at  each  vertex.  If 
a  given  eight-vertex  model  has  a  phase  transition,  then  the  free  energy  per 
vertex  has  a  branch-point  singularity. 

The  eight-vertex  model  has  been  shown  to  be  equivalent  to  two 
Ising  models  with  nearest  neighbor  coupling  interacting  with  another  via  a 
four-spin  coupling  term  (Kadanoff  and  Wegner,  1971:  Wilson  and  Fisher, 

1972).  For  n  >  4,  a  ' 'biconical ' '  fixed  point  describing  a  tetr acri tical 
fixed  point  appears  (Nelson,  Kisterlitz  and  Fisher,  1974) .  An  extension  to 
the  Weiss  model  (Weiss,  1907)  for  ferromagnets  is  also  possible  (Reidel  and 
Wegner,  1969) .  In  this  instance,  the  entropy  S  per  spin  depends  on  the 
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square  of  the  magnetization: 

2  2  ,2 
m  —  mj |  +  m j  , 


to  obtain: 


S  -  s0  -  I  "*  +  •• 


The  free  energy  per  spin: 


F  =  E  -  TS, 


where 


E  =  l/2I0(mf|  +  (1  -  A)m*), 
is  differentiated  to  obtain  the  equations  of  state: 

If  *“||  <*T-I0  +  bT.>)  -k||. 


-  mj^CaT  -  IQ(1  -  A)  +  bTm»)  =  hj^.  (97] 

For  h | |  =0,  the  paramagnetic  region  is  defined  by: 

m|  |  =  0,  (9) 

and  the  ferromagnetic  region  by: 

aT  -  I  +  bTm*  =  0.  (9! 

o 

On  the  A-line  (Seidel  and  Vegher,  1969)  both  conditions  described  by  Equs. 


(98)  and  (99)  hold,  and 

/  hf  \  2 

*x  ■  -  <  r1  ' 

O 


(100) 


where  t  is  the  reduced  temperature  and  or1  =  a/b. 

This  development  may  be  extended  to  include  a  description  of 
ferromagnets  with  a  Hamiltonian  of  cubic  symmetry  (Aharony,  1973),  due  to 
the  presence  in  the  Hamiltonian  of  single  ion  terms  which  reflect  the 
lattice  symmetry,  thus  breaking  the  full  rotational  invariance.  For  d  =  3, 
n  >  3,  a  crossover  from  isotropic  (Heisenberg)  to  characteristic  cubic 
behavior  occurs,  and  there  is  an  approximation  to  the  behavior  of  the  Fez+ 
ion  in  the  hemoglobin  molecule.  The  benefit  to  be  obtained  from  this 


analogy  between  wave  mechanical  systems  and  lattice  systems  is  that  the 
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physics  of  the  latter  predicts  phase  changes  in  the  former  under  appropriate 
conditions. 

The  general  model  may  also  include  the  assumption  that  for 
only  d  -  m  coordinates  is  there  a  contribution  to  the  free  energy  density, 
which  is  quadratic  in  the  spatial  derivatives  of  the  order  parameter,  f 
(Hornreich,  Luban  and  Shtrikman,  1975).  In  this  case,  the  classical-quantum 
mechanical  boundary  is  given  by  d  =  4  +  m/2.  However,  for  clarity's  and 
simplicity's  sake,  we  shall  assume  that  m  =  0  here. 

ii.  Summary  of  ( i) 

The  circumstances  giving  rise  to  the  conditions  described  by 
Equ.  (1)  are  general  and  are  those  of  critical  phenomena. 

The  potential,  <i>  (Equs.  (20),  (63)  and  (76)),  the  minimiza¬ 
tion  of  which  results  in  an  equation  formally  equivalent  to  the  Schrodinger 
equation,  can,  in  another  formal  analogy,  be  identified  with  the  exchange 
energy  at  the  vertex  of  a  lattice.  This  indicates  that  the  one-dimensional 
problem  represented  by  polyacetylene  can  be  generalized  to  n-dimensions . 

The  transition  point  defining  the  boundary  between  single  and 
two  state  behavior  is  then  quite  definite:  the  boundary  is  at  d  =  4.  For 
d  _>  4,  we  have  single  states,  and  for  d  <  4,  we  have  two  states  due  to  the 
breakage  of  symmetry. 

The  principal  quantum  number  and  angular  momentum  quantum 
number  (n)  of  elementary  particle  systems  are  equivalent  to  the 
dimensionality  of  space  (d)  and  the  dimensionality  of  the  spin  space  (n)  for 
n-vector  spins  of  the  general  theory  of  critical  phenomena.  Turning  this 
analogy  around:  we  may  use  the  general  theory  to  indicate  the  transition 
behavior  for  elementary  particle  systems,  ^ad  a  ''Mendeleev's  table1'  may  be 
constructed  based  on  the  d  and  n  parameters.  In  such  a  table,  an  electron 
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system  would  be  mapped  by  the  principal  quantum  number  (d)  and  the  orbital 
angular  momentum  quantum  number  (n).  Such  tables  have  already  been 
constructed  for  lattice  systems  (Fisher,  1974). 
iii.  Application 

The  above  topological  analysis  may  be  extended  to  random 

functions  in  general  of  the  anharmonic  oscillator,  indicating  the  relation  of 

fluctuation  phenomena  ,.0  wave  events  in  general. 

In  the  sections  above,  f  corresponds  with  the  momentum,  p, 

2 

and  thus  for  the  harmonic  oscillator,  f  corresponds  to  coordinate  x  , 

2 

because  E  =  hf  for  quantal  signals,  and  E  =  l/2kx  in  the  case  of  the 
harmonic  oscillator. 

Commencing  with  the  observation  that  when  a  particle  of  mass 

3 

m  and  coordinate  x  acts  under  the  influence  of  force  F(x)  =  -  x  -  ax  -  b, 

the  equation  of  motion  for  such  an  object  is: 

mi  -  qx  =  -x^  -  ax  -  b,  (101) 

where  q  is  a  damping  constant,  and  a  and  b  are  proportionality  constants. 

For  convenience,  we  define  a  unit  mass,  i.e.,  m  -  1,  and  denote  a/m  and  b/m 

by  aQ  and  bQ,  respectively. 

The  potential  energy  of  the  particle  is  given  by: 

V(x)  =  1/ 4x4  +  l/2a  x*  +  b  x,  (102) 

o  o 

as  F(x)  =  -3 V/3x.  Considering  now  the  case  of  local  stability,  we  set  x  =  x 
=  0,  or 

x3+ax+b  =0.  (103 ) 

o  o 

The  positions  of  stability,  which  are  determined  by  the  parameters 
a  and  b  ,  are  thus  obtained  from  Equ.  (103). 
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Fig.  30  Potential  V(i)  =  l/4x  +  l/2ax  +  bx  for  b  =  0  and  a  <  0 
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Mathematically,  a  stable  position  satisfies  the  condition  of 

minimum  potential  energy.  For  the  stable  positions  x^  and  x^  and  b  =  0, 

that  minimum  would  be  V(x.  .)  =  l/4a  a,  and  in  the  case  of  a  particle 

1 1  z  o 

confined  to  the  area  of  instability,  3Va/3xa  >  0,  and  the  particle's 

position  would  be  defined  by  x  > ±/  -aQ/3 . 

We  assume  now  that  the  particle  is  located  near  the  origin 

and  its  fluctuations  into  either  position,  x^  and  x^,  describable  by  a 

random  impulse  train: 

00 

$(t)  =  l  (-1)J8  6<t  -  t'),  (104) 

i=l  i 

where  are  the  magnitudes,  J  is  either  zero  or  one,  and  5(t  -  t*)  is  the 
Dirac  delta  function.  The  correlation  function  of  the  impulse  is  then: 

<£(t)  5(t’>  =  cQ8(t  -  t'),  (105) 

where  cq  is  the  magnitude  of  the  impulse  at  (t  -  t'),  the  mean  time  between 
impulses.  It  follows  that  the  probability  of  finding  the  particle  between 
coordinate  x  and  x  +  dx  at  time  t  is  described  by  the  stochastic  equation 
(Sveshnikov,  1968): 

OO 

*  y  I  i  j 

■f-  (X)  =  j!(-3x)  [kj(x)f(x)].  (106) 

Since  we  are  interested  here  in  the  simple  case  of  continuous 
Markov  processes,  we  set  k^  =  k^  =  ...  =0,  and  obtain  from  Equ.  (106)  the 
Fokker-Planck  equation: 

i  21 

^(x)  =  -3x  [k^x^x)  +  1/2  2  [k2(x)^(x)],  (107) 

dx 

where  k^(x),  k^(x)  are  the  drift  coefficient  and  the  diffusion  coefficient 
respectively.  In  the  present  case,  k^(x)  is  the  force  defined: 
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k  (x)  *  -x*  -  a  x  -  b  =  -  (108) 

1  o  o  ox 

and  the  magnitude  of  the  impulse,  >  is: 

k-  (x)  -  c  .  (109) 

2  O 

We  note  in  passing  that  this  Fokker-Planck  equation  has  been 
related  to  birth-and-death  decriptions  of  fluctuations  (Nicolis  and 
Prigogine,  1977:  Nitzan,  Ortoleva,  Deutch  nd  Ross,  1973). 

We  now  impose  the  assumption  that  the  distribution  function 
is  normalizable  and  vanishes  at  infinity.  With  the  particle  in  the 
stationary  state,  -f-(x)  =  0,  and  substituting  Equs .  (108)  and  (109)  into 
(107),  we  may  solve  for  S  with  respect  to  this  stationary  state  and  obtain: 

-f-U)  -  Nexp  ( 1/ 4x4  +  l/2anx»  +  bnx)},  (110) 

co  0  0 

where  N  is  a  normalization  constant. 

In  the  case  of  b^  *  0,  the  distribution  function.  S',  can  be 

expressed  by  a  Gaussian  wave  form: 

*X)  ■  N0.„  (  -  (HD 

where  nn  =  Nexp  (-0- )  and  a_  and  cft  are  mean  and  variance,  respectively, 

u  zc  u  u 

Equ.  (Ill)  may  now  be  compared  with  the  previous  result  (Equs.  (69))  if  the 
following  replacements  are  made: 

xa  -*■  f,  a^  /a  =  fg  and  1/2Cq  -*■  c(n)2. 

Transferring  the  distribution  function,  S(x) ,  into  the 
frequency  domain  gives: 

S(f)  =  nA  exp{-c(n)2(f  -  f  )*},  (112) 

0  o 

previously  obtained  (Equ.  (71)). 

Considering  now  a  particle  fluctuating  from  stable 
coordinates  xg  to  a  new  coordinate  x^  +  5^,  then  if  the  higher  orders  of  8x 
are  discarded,  a  new  distribution  function  at  x$  +  6x  may  be  calculated  from 
Equ.  (110): 
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f(x  +  Sx)  =  N  exp[-{  —  (l/4x4  +  l/2a  x*  +  b  x  )  (113) 

s  o  '  c  s  osos 

o 

+  7-  (3x*  +  a  >  3c)*}]  . 
c  o 

o 

The  particle  is  then  confined  to  a  new  position  of 
instability  described  by  the  second  term  of  the  exponent  in  Equ.  (13): 

-(3x*  +  aQ)/co.  On  the  other  hand,  from  the  V(x),  x  relation,  we  know  that 

aay  d*v 

=  3x|  +  aQ.  By  imposing  the  condition  of  instability,^—  =  x*  +  aQx  + 

X=X  x=x 

s  s 

bQ  =  0,  we  arrive  at  the  equation  for  the  cusp  in  the  a^,  bQ  plane:  4a*  + 

27b*  =  0,  as  shown  in  Fig.  31.  The  particle  may  fluctuate  on  either  the 

positive  or  negative  side  of  the  x  axis  and  along  the  a  axis  to  a  =0.  If 

o  o 

the  restriction  b^  =  0  is  lifted,  the  particle  may  fluctuate  on  either  side 

of  the  b  axis  and  along  the  a  axis  to  a  =0. 
o  00 

The  circumstance  of  isomerization  in  transpoly acetylene  (or 
isomerization  in  general)  is  thus  related  to  symmetry  breaking  in  critical 
phenomena.  This  is  a  necessary,  but  insufficient,  condition  for  dynamic 
propagation  to  occur.  A  further  condition  is  that  an  energy  gap  be  created 
between  the  isomerized  states.  The  creation  of  this  gap  was  addressed  by 
Hubbard  (1963,  1964a,b). 


Fig.  31  Fluctuation  limits  of  particle:  in  the  x,  b  plane  (x*  +  a  x  +  b 

a  <  0:  b  <  +  1.31):  in  the  x,  a  plane  (-a  +  3x*):  and  in  the  a,  b°plane° 

($a*+27g»).  0  0 

o  o 

XI  THE  HUBBARD  MODEL 

In  sections  I  -  IX  soliton  propagation  was  related  to  exci ton-phonon 
coupling  in  a  one -dimensional  system  with  incommen surat ion.  The 
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incommen surat ion  permitted  an  energy  gap  to  develop  which  established  the 
conditions  for  parametric  excitation  diffusion.  The  soliton  propagation  was 
treated  microscopically,  i.e.,  each  link  in  the  one-dimensional  chain  was 
allocated  a  pump,  signal  and  idler. 

Although  this  approach  is  necessary  and  provides  physical  insight,  a 
macroscopic  picture  provides  empirically  a  simpler  treatment.  In  essence, 
this  means  the  replacement  of  a  many-body  system  by  a  system  of  independent 
quasiparticles  moving  in  a  properly  determined  field  which  may  be  random 
(Economou,  Vhite  and  DeMarco,  1978).  Just  as  Hubbard  proposed  substituting 
for  a  partly  filled  d-band  of  noninteracting  electrons  the  total  spin  of  the 
atom  based  on  correlated  electron  motions,  so  one  may  substitute  for  the 
individual  monomer  units  of  a  one— dimensional  conductor,  a  macroscopic 
picture  of  the  whole  chain,  or  at  least  chain  parts  -  if  correlations  of 
exci ton-phonon  coupling  are  strong  enough.  Although  one  may  still  suppose 
the  microscopic  soliton  to  move  from  monomer  to  monomer  as  described  in 
section  I  —  IX,  its  motion  may  be  correlated  in  such  a  manner  as  to  give 
properties  of  a  macroscopic  theory.  Such  correlations  are  describable  in 
terms  of  Zubarev's  (1960)  Green's  function  technique. 

Hubbard's  model  may  be  viewed  quite  generally.  Just  as  Hubbard  had  in 
mind  the  case  of  d— electrons,  yet  addressed  the  case  of  an  s— band  having  two 
states  per  atom,  so  we  shall  keep  in  mind  the  two  isomerized  states  of 
transpolyacetylene,  while  yet  addressing  the  general  case  described  by 
Equ.  (62b)  which  is  Equ.  (20). 

The  system  addressed  by  Hubbard  behaves  as  though  it  has  two  energy 
levels  Tq  and  T^  +  containing  1  -  l/2n  and  l/2n  states  per  atom,  respec¬ 
tively,  where  I  is  the  interaction  energy.  Here,  on  the  other  hand,  a 
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system  of  two  states  are  addressed,  both  of  energy  Tq,  but  for  which  energy 
I  is  needed,  if  it  is  to  pass  from  one  state  to  the  other.  This  I,  then,  is 
the  applied  current,  I,  of  Table  I  of  section  II  or  the  pumping  energy. 

Given  these  identifications,  there  is  a  direct  analogy  with  the  s-band 
problem  and  the  results  of  the  previous  analysis  by  Hubbard  may  be  adopted. 

Proceeding  directly,  then,  to  certain  conclusions  drawn  from  the 
Hubbard  formulation:  if  P  is  the  density  of  states  of  the  band  structure 
(e.g.,  of  the  atom  or  of  the  whole  (CH)^  chain),  and  p  is  the  density  of 
states  of  the  microsystems  (e.g.,  the  electrons  or  the  (CH)^  monomers),  then 
the  function  g(E)  which  transforms  the  microdensity  into  the  macrodensity, 
p(E)  -  P(g(E) ) ,  is  (Hubbard,  1963,  p.  252): 


PROFILE  OF/>(E) 


Figure  32:  Typical  g(E)  curve.  From  Hubbard,  1964,  p.  253. 
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By  this  relation  we  see  that  the  band  is  split  into  two  parts  by  the  infinity 

of  g(E)  at  E  =  T  +  1(1  -  n  ),  where  ?  =  ct  c .  and  ct  and  c.  are  creation 
o  -<y  ic  to  in  10  i<r 

and  destruction  operators  for  electrons  -  or  conformational  states  in  the  case 

of  (CH)^  -  in  the  state  (i,o)  a  =  +  1. 

As  in  the  limit,  I  -*■  0,  g(E)  is  the  straight  line  A0B,  the  splitting  of 

the  microdensity  of  states  is  due  to  the  perturbation  I.  This  perturbation 

thus  performs  the  same  function  as  the  Peierls  distortion  in  the  (CH)^  system 

or  the  torsion-bar  spring  constant,  K,  of  Table  I  in  section  II. 

Some  correspondences  with  the  Hubbard  formulation  may  now  be  indicated.  If 

<j)^(x>  are  the  degenerate  wave  functions  of  (CH)^,  p  is  a  label  distinguishing 

various  degenerate  wave  functions,  c.^  is  the  destructive  operator  for  a 

soliton  in  the  state  p  on  the  monomer  labeled  i,  i.e.,  for  a  soliton  in  the 

state  <{>  (x-/R.)  where  R.  is  the  position  vector  of  monomer  i  and  cT  is  the 
P  ~  “ i  —  i  ip 

corresponding  creation  operator,  then  for  narrow  energy  bands  the  Hamiltonian 
is  (Hubbard,  1964,  p.  240): 


H  = 


pv 


ij  in  jv 


+  1/2 


i,  j,k,£p,  V<J,T 


(  i,  P,  j  l1/rlka'jlT)ciflcjVc2,TctCT'  (114) 


where 

T^  =  /  <|>*(x-  R.  )t^  V*  +  V]*, ,(x-  R.  )dx, 

ij  p  -i  2m  v  -j 

and 


(115) 


(ip,  jv|l/r|ko,JlT)  =  e1/<ti*(x-R^)({)q(x-R^)(j)^(x1-R  (x^R^jdxdx1  (116) 


X  -  x 


Recalling  Equ.  (113)  above,  noting  the  formal  similarity  of  T^  (Equ. 

(115))  and  N  exp[-{“  (l/4x4  +  l/2a  x*  +  b  x  )] ,  the  second  term  in  Equ.  (114) 
orc0s  osos  H 

is  equated  to  ~(3x|  +  a  Mix)1.  The  general  result  of  section  X  (iii)  above 
cQ  o 

is  then  applicable  to  the  Hubbard  model. 
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It  is  instructive  to  dissect  the  dynamics  of  the  Hubbard  model  further. 
For  example,  if  we  equate  which  are  the  band  energies  of  the  Hubbard 
electron  model  with  the  conformational  energies  of  the  (CH)^  chain,  E^.  with 
the  conformational  energies  of  the  separate  (CH)^  monomers  and  equate: 

(117) 


p, q  E+E  -E  pq  '  TJ+' 
p  q 

where  X  =|lp)(q|  ,  with  the  transition  probabilities  from  one  to  the  other  of 

pq 

CC  1 

the  isomerized  states  of  (CH)  ,  p  and  q,  then  the  function  F(E)  =  e, -T  = 

x  k  o  D(E) 

provides  a  relation  of  the  microsystems  (monomers)  to  the  macrosystem  (the 
chain).  This  relation  exists  between  the  density  of  states  P(E)  of  the 
macrosystem  and  the  density  of  states  p(E)  of  the  microsystems  (Hubbard, 

1964a,  p.  250):  ATOMIC  RESONANCES 


Figure  33:  The  function  F(E)  with  relation  of  the  density  of  states  p(E)  and 
P(E) .  (Hubbard,  1964a,  p.  250). 


Clearly,  then  F(E)  is  the  x  of  Figures  17  and  22A  of  section  IX  and  the  x  of 
Figure  31  above.  As  stated  in  section  X  (iii)  above,  this  x  is  a  substitution 
for  fo  of  Equ.  (62)  and  the  analysis  comes  full  circle. 
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Equation  (62)  is  also  formally  similar  to  the  nonlinear  Schrodinger 
equation  upon  which  the  Davydov  picture  of  vibrational  soliton  propagation 
is  built  (cf.  Davydov  and  Kislukha  (1956)  Equ.  (16):  Davydov  (1980a)  Equ. 
(4.4):  Davydov  (1980b)  Equ.  (2.22):  Davydov  and  Enol'skii  (1980)  Equ. 

(3.2)).  Davydov  subscribes  to  the  notion,  attributed  to  Frohlich,  that  in  a 
one— dimensional  system  the  coupling  of  electrons  with  deformation  may  lead 
to  a  superconductive  effect  without  electron  coupling.  The  possibility  of 
such  coupling  presupposes  the  creation  of  an  energy  gap  and  Figure  32, 
above,  illustrates  the  dynamics  of  the  gap’s  formation. 

In  this  section  missing  pieces  of  the  puzzle  were  provided  concerning 
the  creation  of  that  necessary  gap.  Whereas  soliton  propagation  in  proteins 
may  require  two  (vibrational)  states  separated  by  an  energy  gap,  I,  in  the 
trans  (CH)  system,  on  the  other  hand,  the  two  (conformational)  states  are 
both  of  the  same  energy,  but  separated  by  an  energy  barrier  which  requires 
energy,  I,  to  surmount  (tunnel  through).  Given  this  difference,  the  soliton 
propagation  picture  in  its  dynamics  remains  the  same. 
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XII.  CONCLUSION 

In  section  X  (i)  the  necessary  and  sufficient  conditions  for  the 
obtainment  of  isomerization  in  a  one-dimensional  system  were  described.  In 
section  X  (iii)  the  random  functions  of  an  anharmonic  oscillator  were 
described,  both  for  the  many— body  and  one-body  system.  In  section  XI  the 
one— body  (Hubbard)  system  was  described  and  the  obtainment  of  an  energy  gap 
was  demonstrated  to  be  formally  similar  to  its  obtainment  in  terms  of  the 
anharmonic  oscillator  picture  of  section  X  (iii). 

In  order  to  understand  the  dynamics  of  the  conducting  one-dimensional 
system,  a  parametric  excitation  picture  was  presented  in  sections  I  -  IX. 
The  advantage  of  this  approach  is  a  physical  insight  into  the  dynamics  of 
conduction  in  (CH)^  type  polymers.  The  disadvantage  is  that  a  description 
of  the  whole  chain  in  such  terms  is  intractable.  On  the  other  hand,  in 
sections  X  -  XI  the  previous  many-body  picture  was  related  to  the  one-body 
(Hubbard)  picture.  The  advantage  of  this  is  that  the  whole  chain  becomes 
quite  manageable.  The  disadvantage  is  that  the  dynamics  of  the  problem  is 
only  introduced  in  a  quite  artificial  way  (if  at  all).  Clearly,  therefore, 
both  pictures  are  necessary  for  full  understanding. 

There  is  a  generalization  of  this  application  of  the  double  potential 
well  model  to  conduction  in  general  in  organic  polymers.  The  relationship 
between  one-dimensional  theories  in  field  theory  and  statistical  mechanics 
is  of  importance  in  the  study  of  quasi  one-dimen sonal  systems  due  to  the 
relation  of  the  Bethe  Ansatz  (Bathe,  1931),  on  the  one  hand,  and  the 
transfer  matrix  methods  of  Onsager  (1944)  and  Baxter  (1972)  on  the  other 
(cf.  Fowler,  1981;  Berkan  and  Cooper,  1981).  In  particular,  the 
diagonalization  of  the  Hamiltonian  or  the  transfer  matrix  has  been  achieved 
in  all  these  models  by  one  form  or  another  of  the  Bethe  Ansatz.  Also,  this 
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method  was  used  to  diagonalize  the  massive  Thirring  model  (Bergnoff  and 
Thacker,  1979),  the  Gross-Neveu  model  (Andrei  and  Lowenstein,  1979)  and  the 
Kondo  problem  (Andrei,  1980). 

Furthermore,  recently  Campbell  and  Bishop  (1981,  1982a, b;  Bishop,  1980) 
demonstrated  the  relation  of  solitons  in  polyacetylene  to  relativistic- 
field-theory  models  and  established  the  existence  within  a  continuum 
electron— phonon  model,  of  a  polaron— like  excitation  additional  to  and 
different  from  the  amplitude  ''kink''. 

The  relativistic— field-theory  aforementioned  is  the  Gross-Neveu  model 
(Gross  and  Neveu,  1974)  described  by  a  Lagrangian  density  of  anharmonic 
oscillator  type.  Polyacetylene  exemplifies  the  situation  of  a  single 
electron  spectrum  with  two  states  symmetrically  placed  in  the  gap  at  c—  = 

■“*  »  which  exists  only  when  there  is  a  single  unpaired  electron  occupying 
the  £  state,  or,  when  the  £  state  is  only  singly  occupied  by  a  localized 
hole  state.  Since  the  electron  or  hole  is  trapped  in  a  structure  due  to 
phonon  coupling,  the  excitation  is  polaron-like.  Similar  results  have  been 
obtained  by  Brazovskii  (1978,  1980)  and  Takayama,  Lin-Liu  and  Haki  (1980). 
This  continuum  model  has  also  been  addressed  by  Horovitz  (1980,  1981)  and 
Krumhansl,  Horovitz  and  Heeger  (1980). 

These  developments  differ  from  the  Su,  Schrieffer  and  Heeger  (1979, 
1980)  treatment  of  electrical  conduction  in  polyacetylene  in  that  the 
soliton— like  polarons  are  obtained  from  continuum  models,  whereas  Su  et  al., 
imposed  a  hyperbolic  tangent  (tanh)  kink  form  on  the  Pierls  amplitude  order 
parameter  to  obtain  kink  solutions.  Both  Brazovskii  and  Takayama  et  al.,  on 
the  other  hand,  demonstrate  that  the  order  parameter  profile  is  an  exact 
solution  in  the  continuum  limit  of  the  self-consistent  equations  of  motion 
for  phonon  and  electron  operators  determined  in  a  variational  scheme.  Thus, 
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this  development  is  fully  compatible  with  the  physical  scheme  represented 
here. 

Campbell  and  Bishop  also  point  out  an  exact  formal  equivalence  between 
the  Brazovskii— Takayama  et  al.,  continuum  equations  and  the  static  classical 
equations  of  the  N=2  Gross-Neveu  model  (Gross  and  Neveu,  1974;  Dashen, 
Hasslacher  and  Neveu,  1974c;  1975a, b).  By  means  of  the  mean  field 
approximation,  the  Bogliubov— de  Gennes  equations  relating  electron  and 
phonon  fields  were  obtained.  Dashen  et  al.,  demonstrated  that  there  are 
three  classes  of  solutions  to  the  Gross-Neveu  Lagrangian:  (1)  The  ground- 
state  solution;  (2)  The  kink  solution;  and  (3)  The  ''bag''  solution 
involving  the  presence  of  fermions  in  occupied  positive  energy  bound  states 
and  the  absence  of  fermions  from  negative  energy  bound  states  in  the  case  of 
antiparticle  or  hole  excitation.  The  field  trapping  the  fermion  is  driven 
by  phonon  coupling  to  the  fermion,  i.e.,  the  fermion  is  a  polaron. 

Also,  Zamolodchikov  and  Zamolodchikov  (1978,  1919)  demonstrated  how  to 
determine  the  exact  S  matrices  of  two-dimensional  models  including  the  sine— 
Gordon  equations  and  the  (W ) 1  model.  Witten  (1978)  determined  that  the 
Lagrangian  model  used  may  be  considered  an  N— component  Majorana  Fermi  field. 
For  N=2  the  theory  is  of  a  single  interacting  Dirac  field  and  is,  in  fact, 
the  massive  Thirring  model  related  to  soliton  theory  by  Coleman  (1975).  For 
N=3  the  field  is  equivalent  to  a  supersymmetric  sine— Gordon  equation  and  for 
N=»4,  to  two  decoupled  sine-Gordon  systems. 

The  general  anharmonic  oscillator  model  forming  the  basis  of  these 
developments  has  been  most  extensively  studied  by  Dashen,  et  al.  Their 
methods,  in  turn,  are  based  on  those  of  Gurzwiller  (1967,  1970,  1971)  and 
Maslov  (1970).  The  meani.ig  of  the  WEB  approximation  was  examined  from  a 
functional- integral  approach  and  the  conclusion  was  made  that  certain  kinds 
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of  nonlinear  field  theories  have  extended  particle  solutions  which  survive 
quantization.  Theories  with  spontaneous  symmetry  breaking  were  shown  to 
produce  hadron— like  objects.  In  these  instances,  the  field  theory  becomes  a 
model  for  a  superconductor.  Models  in  four  space-time  dimensions  can  be 
constructed  of  SU(2)  Tang— Mills  fields  coupled  to  fermions  and  a 
spontaneous-asymmetry-breaking  scalar  isospinors. 

Dashen  et  al.  (1975a)  also  used  the  WKJhf  method  to  compute  the  particle 
spectrum  of  the  Gross-Neveu  model  and  found  a  spectrum  of  particles  of 
f ermion-ant if ermion  bound  states  and  multifermion  bound  states.  The  close 
similarity  between  the  sine-Gordon  equation  and  the  Gross-Neveu  model  was 
noted,  but  the  reason  for  the  correspondence  not  discerned.  It  is,  however, 
probably  not  unimportant  that  the  Lagrangian  for  both  theories  are  examples 
of  Riemann-Hugoniot  singularities.  Dashen  et  al  (1975b)  pointed  out  the 
correspondence  of  the  sine-Gordon  coupling  parameter  and  the  four  fermion 
coupling  of  the  Thirring  model.  The  fermions  are  thus  solitons. 

In  conclusion,  therefore,  the  enharmonic  oscillator  model,  whether  in 
the  form  of  linked  pendula  model,  or  conducting  polymer  form,  offers  a  rich 
spectrum  of  particles.  Due  to  the  so-called  ''democracy''  of  the  sine- 
Gordon  equation  (Dashen  et  al  1975b),  there  is  no  fundamental  distinction 
between  the  particle  of  the  n  =  1  quantum  state  and  the  higher  mass  n  >  1 
particles.  As  the  anharmonic  oscillator  model  indicates  electron-phonon 
coupling  involvement  in  electrical  conduction,  and  the  electron  can  occupy 
any  of  those  states,  one  may  expect  a  rich  variety  of  polymers  able  to 
conduct  either — electrically  or  vibrat ionally — besides  the  one-dimensional. 
This  is  because  application  of  the  anharmonic  oscillator  model  is 
appropriate  under  any  condition  in  which  the  double  potential  well  model 
applies. 
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